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EXISTENCE, CHARACTERIZATION AND STABILITY 
OF PANSU SPHERES IN SUB-RIEMANNIAN 3-SPACE FORMS 

ANA HURTADO AND CESAR ROSALES 


Abstract. Let M be a complete Sasakian sub-Riemannian 3-manifold of constant Webster scalar 
curvature re. For any point p £ M and any number A £ R with A 2 + re > 0, we show existence 
of a C 2 spherical surface S\(p) immersed in M with constant mean curvature A. Our construc¬ 
tion recovers in particular the description of Pansu spheres in the first Heisenberg group [48] and 
the sub-Riemannian 3-sphere [35]. Then, we study variational properties of S\(p) related to the 
area functional. First, we obtain uniqueness results for the spheres S\(p) as critical points of the 
area under a volume constraint, thus providing sub-Riemannian counterparts to the theorems of 
Hopf and Alexandrov for CMC surfaces in Riemannian 3-space forms. Second, we derive a second 
variation formula for admissible deformations possibly moving the singular set, and prove that 
S\ (p) is a second order minimum of the area for those preserving volume. We finally give some 
applications of our results to the isoperimetric problem in sub-Riemannian 3-space forms. 


1. Introduction 

The study of variational questions related to the area functional in sub-Riemannian geometry is 
a focus of recent attention, especially in the Heisenberg group HP, which is the simplest non-trivial 
example of a flat sub-Riemannian manifold. Indeed, the isoperimetric problem in HI", where we seek 
regions minimizing the (sub-Riemannian) perimeter under a volume constraint, is one of the main 
open problems in this field with several important contributions in the last ten years. 

The existence and boundedness of isoperimetric solutions in H" are consequence of more gen¬ 
eral results of Leonardi and Rigot for Carnot groups [39]. Moreover, the first variational formulas 
for area and volume imply that the boundary of any C 2 solution is a compact hypersurface with 
constant mean curvature (CMC) in sub-Riemannian sense, see [13] and [53]. In 1982 it was con¬ 
jectured by Pansu [48] that the isoperimetric regions in H 1 are congruent to the topological balls 
enclosed by a family {Sa}a>o °f CMC spheres of revolution with C 2 regularity, see [45], [38] and 
[53] for explicit expressions of S\. These spheres are geometrically obtained by rotating about the 
vertical axis in H 1 a Carnot-Caratheodory geodesic (CC-geodesic) connecting two points on this 
axis. Such a geodesic is a Euclidean helix with curvature depending only on the vertical distance 
between its extreme points. Pansu’s conjecture has been supported by many partial results, where 
further hypotheses involving regularity and/or symmetry of the solutions are assumed. We refer the 
reader to [8, Ch. 8] and the Introduction of [52] for a precise account and description of the related 
works. More recently, Monti [46] has analyzed symmetrization in B n , and Cheng, Chiu, Hwang and 
Yang [10] has used a notion of umbilicity to characterize Pansu’s spheres in HP with n ^ 2. It is 
worth mentioning that the optimal regularity for isoperimetric sets in HP is still an open question. 
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Indeed, though the conjectured solutions are C 2 smooth, there exist non-compact area-minimizing 
surfaces in H 1 with less regularity, see [49], [15], [47] and [51]. 

Besides the Heisenberg group, it is also natural to investigate the isoperimetric question in other 
sub-Riemannian 3-spaces, especially in those representing the most simple and symmetric exam¬ 
ples of the theory. These are the sub-Riemannian 3-space forms , defined in Section 2.2 as com¬ 
plete Sasakian sub-Riemannian 3-manifolds of constant Webster scalar curvature. By a result of 
Tanno [60] there is, up to isometries, a unique simply connected sub-Riemannian 3-space form M(/t) 
of Webster scalar curvature «. The space M(k) is the Heisenberg group H 1 for k = 0, the group 
of unit quaternions § 3 C R 4 for k = 1, and the universal cover of the special linear group SL(2,R) 
for k = —1. It is also known that the spaces M(re) are the unique simply connected, homogeneous, 
contact sub-Riemannian 3-manifolds with the largest possible isometry group, see for instance [21]. 
For an arbitrary sub-Riemannian 3-space form M, standard Riemannian arguments show that M is 
isometric to M(k)/G, where k is the Webster scalar curvature of M and G is a subgroup of isometries 
of M(k) isomorphic to the fundamental group of M. This permits the construction of 3-space forms 
with non-trivial topology like flat and hyperbolic cylinders I 2 x S 1 , as well as a sub-Riemannian 
model of the real projective space RP 3 , see Example 2.2 for details. 

Our main motivation in this paper is to analyze a possible extension of Pansu’s conjecture in¬ 
side any sub-Riemannian 3-space form M. For that, we must produce first CMC spheres S\{p) 
playing the same role in M as Pansu’s spheres in the Heisenberg group M(0). This was already 
done in the sub-Riemannian 3-sphere M(l), where the authors employed the explicit expression of 
the CC-geodesics to construct S\(p ), see [35]. In the present work we avoid the computation of 
CC-geodesics; instead of that, we will use CC-Jacobi fields (introduced in Section 2.3) systemati¬ 
cally. This approach allows us to introduce and study in a unified way the spheres S\ (p) for general 
sub-Riemannian 3-space forms. 

In Section 3 we define S\(p) as follows. For a fixed point p £ M and a number A £ R we consider 
the flow F[0 , s) := 70 (s) of all the CC-geodesics of curvature A in M leaving from p with unit hori¬ 
zontal velocity 6. Thus, in order to produce a spherical surface, we require the CC-geodesics jg to 
meet again in M. This is clearly equivalent to that the associated vector fields Vg(s) := ( 8F/d0)(s , 6) 
vanish for some positive value. The remarkable fact is that Vg is a CC-Jacobi field along 70 with 
Vfl(0) = 0 and so, it is explicitly computable since the Webster scalar curvature of M is constant. 
From this calculus we deduce in Lemma 3.1 that the arithmetic condition A 2 + k > 0 is equivalent 
to the existence of a first cut point p\ for the CC-geodesics jg(s) when s := Tt/^/\ 2 + k. This leads 
us to consider the set 


S\(p) := {7e(s); 0 < s ^ 7r/\/A 2 + k}, 

where 9 ranges over the unit horizontal circle of M at p. In the model spaces M(0) and M(l) the 
previous definition recovers the spherical surfaces described in [48] and [35]. In the hyperbolic space 
M(k) with n < 0 we show in Theorem 3.3 that S\(p) is also an embedded C 2 surface congruent to 
a sphere of revolution with C°° regularity off of the poles {p,p\}. So, by taking into account that 
any sub-Riemannian space form M is a quotient of some M(k), we are able to prove in Theorem 3.3 
that S\(p) is a C 2 oriented sphere immersed in M with C°° regularity off of the singular set {p,p\}, 
where the tangent plane equals the horizontal plane. Moreover, in Theorem 3.3 we exhibit some 
geometric properties of S\(p) regarding its symmetry and embeddedness. Precisely, S\(p) is always 
Alexandrov embedded and invariant under the isometries of M fixing p. As a difference with re¬ 
spect to the spaces M(k) the sphere S\{p) need not be embedded in M for small values of A, see 
Remarks 3.5. After some computations in the coordinates (0, s) we deduce that S\(p) has CMC 
equal to A off of the poles. Furthermore, if Nh denotes the horizontal projection of a unit normal 
N along S\(p), then the function |^| _1 turns out to be integrable, see Corollary 3.7. This is a 
non-trivial fact (since |A^| vanishes at the poles) that will be used frequently through this work. 
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Once we have constructed the spheres S\(p) it is natural to analyze if they provide nice candi¬ 
dates to solve the isoperimetric problem in any sub-Riemannian 3-space form M. For that, we first 
observe that S\(p) is volume-preserving area-stationary , i.e., a critical point of the area functional 
for deformations preserving volume. This is an immediate consequence of the characterization re¬ 
sult in [54], see Theorem 4.5, since S\(p) is a C 2 oriented CMC surface with finitely many singular 
points. Now, an interesting question is to characterize S\ (p) in the family of volume-preserving area 
stationary surfaces, possibly under additional hypotheses. In Section 4.2 we prove some uniqueness 
results in this line providing sub-Riemannian counterparts to the well-known Hopf and Alexandrov 
theorems for CMC surfaces in Riemannian manifolds. 

Recall that Hopf’s theorem [32] establishes that a topological sphere immersed in R 3 with CMC 
is a round sphere. This was extended by Chern [17] to arbitrary complete Riemannian 3-manifolds 
of constant sectional curvature. In Theorem 4.7 we show that a C 2 oriented sphere with CMC 
immersed in a sub-Riemannian 3-space form M must be a spherical surface S\(p). Recently Cheng, 
Chiu, Hwang and Yang [10] have obtained a Hopf theorem for umbilic hypersurfaces in the Heisen¬ 
berg group H” with n ^ 2. On the other hand, Alexandrov’s theorem [1] states that a compact CMC 
surface embedded in R 3 is a round sphere. The proof, based on the so-called Alexandrov reflection 
technique, can be readily generalize to the hyperbolic space and a hemisphere of § 3 . In Theorem 4.11 
we characterize the spheres S\{p) as the unique compact volume-preserving area-stationary surfaces 
of class C 2 in M(re) with k ^ 0, or inside an open hemisphere of M(k) with k > 0. This theorem was 
previously established in [54] for the particular case of M(0). The vertical Clifford tori in § 3 illus¬ 
trate that the Alexandrov theorem is no longer true if the surface is not entirely contained inside an 
open hemisphere. As we see in Remark 4.12 the theorem also fails for sub-Riemannian 3-space forms 
with non-trivial topology. The proofs of Theorems 4.7 and 4.11 rely on results in [13] and [54] about 
geometric, topological and variational properties of C 2 critical points of the area under a volume 
constraint, see Section 4.1 for precise statements. These were employed in [54] and [35] to describe 
complete C 2 volume-preserving area-stationary surfaces with non-empty singular set in M(0) and 
M(l). With similar techniques Galli characterized C 2 area-minimizing and area-stationary surfaces 
in the roto-translation group [24] , and in the sub-Riemannian Sol manifold [25] . Recently, Bernstein 
type problems and local minimizers of the area with C 1 or less regularity have been studied in the 
Heisenberg groups H", see the Introduction of [27] for a complete account of related works. 

Once we know that the spheres S\(p) are volume-preserving area-stationary surfaces, we analyze 
in Section 5 if they are also second order minima of the area under a volume constraint , thus pro¬ 
viding natural candidates to solve the isoperimetric problem. As a fundamental tool, we need the 
second derivative of the functional A + 2W stated in Theorem 5.2 for certain deformations of S\(p). 
The second variation of the sub-Riemannian area has appeared in several contexts, see for instance 
[13], [4], [18], [41], [47], [34], [12], [31], [24], [43] and [27]. The derivation of the second variation 
formula involves a technical problem since, in contrast to the first variation, the differentiation under 
the integral sign is not always possible. Nevertheless, this problem is overcome if we require the 
employed deformations to be admissible. These are introduced in Definition 5.1 as those satisfying 
the conditions necessary to apply the usual Leibniz’s rule for differentiating under the integral sign, 
see Lemma 2.10. We must remark that, unlike some previous approaches, our deformations are only 
assumed to be admissible; they are no restricted to fixing the singular set nor having a special form. 
Indeed, the notion of admissibility permits us to obtain a general second variation formula for CMC 
surfaces in arbitrary Sasakian sub-Riemannian 3-manifolds, see Theorem A.l. The proof of Theo¬ 
rem 5.2 comes from Theorem A.l with an additional effort to demonstrate that all the divergence 
terms in the general second variation formula vanish for S\(p). For non-admissible variations the 
second derivative A"(0) may exist or not and, in case of existence, it is not easy to compute it. 
Some particular examples of these situations for CMC surfaces with singular curves were discussed 
in [34], [12] and [24]. 
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From Theorem 5.2 the second derivative of the functional A + 2AU associated to an admissi¬ 
ble variation of a sphere S\ (p) is controlled by a quadratic form X, called the index form, which 
involves the normal component of the velocity vector field and some geometric quantities defined 
off of the poles. Surprisingly, the expression in polar coordinates of I does not depend on the 
ambient manifold M nor the point p £ M, which suggests that all the spheres S\ (p) will share 
the same behaviour with respect to stability conditions. Indeed in Proposition 5.5 we deduce that 
S\(p) is always strongly stable, that is (A + 2AU) ,/ (0) 0, for any admissible variation vanishing 

at the poles or along the equator. This generalizes a previous result of Montefalcone [42] for some 
particular variations of Pansu’s spheres in the Heisenberg group M(0). Nevertheless the spheres 
S\(jp ) are not strongly stable for arbitrary admissible variations. To see this it is enough to consider 
the deformation of S\(p) by Riemannian parallel surfaces, see Remark 5.8 for details. However, 
by using analytical arguments based on double Fourier series, we are able to prove in Theorem 5.9 
that S\ (p) is stable under a volume constraint, i.e., it has non-negative second derivative of the area 
under volume-preserving admissible deformations. In the Heisenberg group H™, Montefalcone [42] 
analyzed a 1-dimensional eigenvalues problem to establish that the Pansu spheres satisfy the same 
stability property for certain radially symmetric deformations. It is worth mentioning that, though 
our stability results are proved for admissible deformations of S\(p), they are still valid for any 
variation for which the second variation formula in Theorem 5.2 holds. Moreover, as it is shown in 
Appendix B, the family of admissible variations of S\(p) is very large. 

We finish the present work with Section 6, where we come back to the isoperimetric problem 
inside any sub-Riemannian 3-space form M. Note that, if M is homogeneous, then we can apply 
the results of Galli and Ritore [26] to ensure existence and boundedness of isoperimetric regions 
for any given volume. Unfortunately, the regularity of the minimizers is still an open question. 
However, it is clear that the boundary of a minimizer must be volume-preserving area-stationary. 
In particular, if S2 is a C 2 isoperimetric region in the model spaces M(k) with k ^ 0, then we can 
invoke Theorem 4.11 to deduce that dfl must be one of the spherical surfaces S\(p). This extends 
to the sub-Riemannian hyperbolic models the uniqueness of C 2 solutions in M(0) proved in [54]. 
In the spherical models M(k) with k > 0 the same conclusion holds for C 2 isoperimetric regions 
contained inside an open hemisphere of § 3 . In [36] we use the second variation formula to discard 
the existence in M(re) with n > 0 of stable tori under a volume constraint, so that, for arbitrary 
C 2 isoperimetric regions, we get the same uniqueness result as in M(k) with k ^ 0. From this we 
may suggest that the spheres S\(p) uniquely minimize the perimeter under a volume constraint in 
M(k). This is a generalization of Pansu’s conjecture in the Heisenberg group M(0). Indeed, our 
stability result in Theorem 5.9 provides further evidence supporting this extended conjecture. On 
the other hand, in a homogeneous space form M of non-trivial topology, a sphere S\(p) need not 
bound a minimizer. This is shown in Example 6.2, where we find vertical tori inside a Sasakian flat 
cylinder 1R 2 x S 1 which are isoperimetrically better than the spheres S\(p) of the same volume. So, 
an interesting question in this setting is the following: are the topological balls enclosed by S\(p) 
isoperimetric solutions for small volumes? 

The paper is organized into six sections and two appendices. The second section contains back¬ 
ground material about sub-Riemannian 3-manifolds. In Section 3 we construct the spherical surfaces 
S\(p) and show some of their geometric properties. In the fourth section we establish uniqueness 
results for S\(p) and describe volume-preserving area-stationary surfaces with non-empty singular 
set. In Section 5 we obtain stability results for admissible variations of S\(p). Some applications 
to the isoperimetric problem in sub-Riemannian 3-space forms are given in Section 6. The second 
variation formula is proved in detail in Appendix A. Finally, in Appendix B we produce several 
examples of admissible variations. 

Acknowledgements. The authors wish to thank the referee for his/her detailed report and useful 
suggestions to improve the presentation of the paper. 
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2. Preliminaries 

In this section we introduce the notation and establish some results that will be used throughout 
the paper. 

2.1. Sasakian sub-Riemannian 3-manifolds. 

A contact sub-Riemannian manifold is a connected manifold M with empty boundary together 
with a Riemannian metric gh defined on an oriented contact distribution R, which we refer to as 
the horizontal distribution. We denote by J the orientation-preserving 90 degree rotation in (R,gn )• 
This is a contact structure on R since J 2 = —Id. 

The normalized form is the contact 1-form g on M such that Ker(?y) = R and the restriction 
of the 2-form dg to R equals the area form on R. We shall always choose the orientation of M 
induced by the 3-form g A dg. The Reeb vector field associated to g is the vector field T transversal 
to R defined by g(T) = 1 and dg(T, U) = 0, for any U. We extend J to the tangent space to M by 
setting J(T) := 0. A vector field U is horizontal if it coincides with its projection Uh onto R. If U 
is always proportional to T then we say that U is vertical. 

The canonical extension of gr is the Riemannian metric g = (•, •) on M for which T is a unit 
vector field orthogonal to R. For this metric it is easy to check that 

(2.1) (J(U),V) + (U,J(V)) = 0, 

for any U and V. We say that M is complete if (M, g) is a complete Riemannian manifold. 

By a local isometry between contact sub-Riemannian 3-manifolds we mean a smooth ( C °°) map 
(j> : M —> M' whose differential at any p £ M is an orientation-preserving linear isometry from R p 
to R'^ p y These maps preserve the normalized forms, the Reeb vector fields, the associated com¬ 
plex structures and the canonical extensions of the sub-Riemannian metrics. An isometry is a local 
isometry which is also a diffeomorphism. A contact sub-Riemannian 3-manifold M is homogeneous 
if the group Iso (M) of isometries of M acts transitively on M. Such a manifold must be complete. 
We say that two subsets Si and S 2 of M are congruent if there is (j) £ Iso(M) such that (f>(Si) = S' 2 . 

By a Sasakian sub-Riemannian 3-manifold we mean a contact sub-Riemannian 3-manifold M 
where gh is a Sasakian metric, this means that any diffeomorphism of the one-parameter group as¬ 
sociated to T is an isometry. This implies, see [6, p. 67, Cor. 6.5 and Thm. 6.3], that the Levi-Civita 
connection D associated to g satisfies the equalities 

D V T = J(U), 

(2.2) D v (J(V)) = J(DuV) + (V, T) U - (U, V) T, 

for any vector fields U and V on M. In particular, any integral curve of T is a geodesic in (M, g) 
parameterized by arc-length. We refer to these curves as vertical geodesics. 

Now we recall some curvature properties of a Sasakian 3-manifold. The curvature tensor in (M, g) 
is defined by 

R(U, V)W := D v DuW - DuD v W + D [uy] W ,, 
where [U,V] is the Lie bracket. From [6, Prop. 7.3] and the fact that (U,T) = g(U), we have 

(2.3) R(U, V)T = (U,T) V — (V,T) U. 

The Webster scalar curvature K of M is the sectional curvature of R with respect to the Tanaka 
connection [6, Sect. 10.4], If Kh denotes the sectional curvature of R in ( M,g ), then 

K=\{K h + 3), 
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which is clearly invariant under local isometries. The following formulas valid for U horizontal and 
V arbitrary were proved in [55, Lem. 2.1] 

(2.4) R(U, V)U = (4 K - 3) < V ., J{U)) J{U) + \U\ 2 (V, T) T , 

(2.5) Ric(V, V) = (4 K - 2 ) |Vfc | 2 + 2 < V,T)\ 

where Ric is the Ricci curvature in (M, g) defined, for any two vectors fields U and V, as the trace 
of the map W H- R{U , W)V. 

Let M be a contact sub-Ricmannian 3-manifold. The homothetic deformation of gh is the one- 
parameter family of sub-Riemannian metrics {{gh)e}s>o on di defined by 

( 2 . 6 ) (g h ) e := e 2 g h . 

This induces new structures of contact sub-Riemannian 3-manifold on M, see [ 6 , p. 103] and [59]. 
The normalized form rj E , the Reeb vector field T e , the complex structure J £ and the canonical 
extension g e are, respectively, given by r) E = e 2 77 , T e = T/e , J E = J , and 

(2.7) g E (U, V) = e 2 g(U , V) + e 2 {e 2 - 1) g(U, T) g(V, T). 

It is easy to check that any isometry for gh is also an isometry for (gh)e- Therefore ( gh)e is a 
Sasakian metric if and only if the same holds for gh- Moreover, in such a case, we obtain 

I< 

K e = 

E z 

where K e is the Webster scalar curvature for ( gh) E . 

2.2. Three-dimensional space forms. 

For k = —1,0,1, we denote by N(re) the complete, simply connected, Riemannian surface of 
constant sectional curvature 4k described as follows. If k = 1, then N(l) is the unit sphere § 2 with 
its standard Riemannian metric scaled by 1/4. If k = —1,0, then N(k) := {p£i 2 ; \p\ < 1/|k|} 
endowed with the Riemannian metric p 2 (dx 2 + dy 2 ), where p(x : y) := (1 + n(x 2 + y 2 )) -1 . Hence 
N(—1) is the Poincare model of the hyperbolic plane and N(0) is the Euclidean plane. 

For k = —1,0 we denote M(k) := N(k) xR. If k = 1, then we set M(l) = § 3 . The space M(k) 
admits a structure of complete Sasakian sub-Riemannian manifold of Webster scalar curvature k, see 
[55, Sect. 2.2] and [ 6 , Sect. 7.4] for an explicit description. The sub-Riemannian metric gh in M(k) is 
of bundle type in the sense of [44, p. 18]. In fact, there is a Riemannian submersion T : M(k) —X N(k) 
for which R = (Kei(dJ : )) ± . Here the orthogonal complement is taken with respect to the canon¬ 
ical extension g of gh. The map J- is the Euclidean projection F(x,y,t) := (x,y) if k = —1,0, 
and the Hopf fibration J 7 (xi,yi,x 2 ,y 2 ) ■= (*? + y\ ~ x\ - yf, 2 {x 2 yi - Xiy 2 ), 2 ( 2 : 12:2 + yiy 2 )) if 
k = 1. We shall always consider the orientation on N(k) for which the restriction of dT to R is 
an orientation-preserving isometry. The oriented fiber of J- through p is the vertical geodesic in 
(M(k),p) parameterized by p + tT p if k = —1,0 or e lt p if k = 1. We call vertical axis the fiber 
through o := (0,0,0) if k ^ 0 or through o := (1,0,0, 0) if k > 0. The diffeomorphisms associated 
to T are translations along the fibers. By means of the homothetic deformation in (2.6) we can 
construct, for arbitrary k € R, the Sasakian space M(k) of constant Webster scalar curvature k. 

In M(k) there is a product * such that (M(k), *) is a Lie group. More precisely, we get the Heisen¬ 
berg group when k = 0 , the group of unit quaternions when k > 0 , and the universal covering of the 
special linear group SL(2,M) when k < 0. The associated left translations (resp. right translations) 
are isometries of M(k) when k ^ 0 (resp. k > 0). Hence the space M(k) is homogeneous and, in 
particular, complete. Any Euclidean vertical rotation (a rotation about the vertical axis) is also an 
isometry of M(k). 

By a result of Tanno [60] the space M(k) is, up to isometries, the unique complete, simply con¬ 
nected, Sasakian sub-Riemannian 3-manifold with constant Webster scalar curvature k. This result 


PANSU SPHERES IN SUB-RIEMANNIAN 3-SPACE FORMS 


7 


together with a standard argument involving Riemannian coverings (see the proof of Prop. 4.3 in 
[20, Chap. 8 ] for details) allows to deduce the following fact. 

Proposition 2.1. If M is a complete Sasakian sub-Riemannian 3-manifold of constant Webster 
scalar curvature k, then there is a surjective local isometry II : M(k) — > M. Moreover, there is 
a subgroup G of isometries of M(k) such that M is isometric to the quotient M(k)/G, and the 
fundamental group of M is isomorphic to G. 

In the sequel, by a 3-dimensional space form we mean a complete Sasakian sub-Riemannian 
3-manifold of constant Webster scalar curvature. 

Example 2.2. The previous proposition suggests that one may construct 3-space forms with non¬ 
trivial topology by taking a non-trivial subgroup G of isometries of M(k) such that the Sasakian 
structure and the group operation in M(re) descend to the quotient M := M (k)/G. In M(k) with 
k ^ 0 this can be done if we consider the subgroup G = Z of Iso(M(k)) given by G := {<p n ; n G Z}, 
where (f> n (p) '■= P + (2nir) T p . The resulting space form is homeomorphic to the cylinder R 2 x S 1 . 
This construction also holds in M(k) with n > 0 for the subgroup G = Z 2 of Iso(M(k)) given 
by the identity map and the antipodal map in § 3 . This leads to a sub-Riemannian model of the 
3-dimensional real projective space RP 3 . 

For a piecewise smooth curve a : I —> N(k) we define a lift of a as a piecewise smooth curve 
5 : I M(k) such that T o 5 = a. From basic properties of principal bundles [37, p. 88 ] there is 
a unique lift with horizontal tangent vector wherever it exists ( horizontal lift) and passing through 
a given point of M(k). The next lemma states that the holonomy displacement associated to a 
horizontal lift of a Jordan curve in N(k) coincides with the area enclosed by the curve. The cases 
k = 0,1 are found in [44, Sect. 1.3] and [50, Proof of Prop. 1]. In general, the proof follows by 
adapting the arguments in [58, Thm. 1, p. 191]. 

Lemma 2.3. Let D C N(re) be a disk of area A bounded by a piecewise smooth Jordan curve. Let 
a : [0,L] — > N(k) be a parameterization by arc-length of dD which reverses the orientation induced 
by D. If a : [0,L] —> M(k) is a horizontal lift of a, then the oriented distance between 5(0) and 
a(L) along the fiber of T through 5(0) equals 2A. 

2.3. Carnot-Caratheodory geodesics and Jacobi fields. 

Let M be a Sasakian sub-Riemannian 3-manifold. A horizontal curve in M is a C 1 curve 7 whose 
velocity vector 7 lies in “H. The length of 7 over a compact interval / C R is given by J T | 7 (s)| ds, 
where | 7 (s)| denotes the sub-Riemannian length of j(s). The Carnot-Caratheodory distance between 
p, q G M is the infimum of the lengths of all C°° horizontal curves joining p and q. Since R is a 
bracket generating distribution such curves exist by Chow’s connectivity theorem [28, Sect. 1.2.B]. 

Suppose that 7 is a C 2 horizontal curve parameterized by arc-length. We say that 7 is a Carnot- 
Caratheodory geodesic, or simply a CC-geodesic, if 7 is a critical point of length under C 2 variations 
by horizontal curves. Following variational arguments, see for example [54, Prop. 3.1], we deduce 
that 7 is a CC-geodesic if and only if there is a constant A 6 R, called the curvature of 7 , such that 
7 satisfies the second order ODE 

(2.8) 7 ' + 2AJ(7) = 0, 

where the prime ' stands for the covariant derivative along 7 in (M, g). If p £ M and v € R p with 
|u| = 1, then the unique solution 7 to (2.8) with 7 ( 0 ) = p and 7 ( 0 ) = v is a CC-geodesic of curvature 
A since the functions ( 7 , T) and | 7| 2 are constant along 7 . The fact that any CC-geodesic is defined 
on R is equivalent to the completeness of M by [5, Thm. 1.2]. 

The next result shows the behaviour of CC-geodesics with respect to local isometries and the 
homothetic deformation in (2.6). The proof comes from (2.7) and the Koszul formula [20, p. 55]. 
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Lemma 2.4. Let M be a Sasakian sub-Riemannian 3-manifold, and 7 : / —> M a CC-geodesic of 
curvature A defined over an open interval I CM. Then we have 

(i) If 4> : M -A M' is a local isometry then </> o 7 is a CC-geodesic of curvature A in M'. 

(ii) For any e > 0, the curve 7 e : I e —> M defined on I E := el by 7 E (s) := 7(s/e) is a CC-geodesic 

of curvature A/e for the sub-Riemannian metric {gh)e ~ £ 2 9h- 

In a model space M(k) any CC-geodesic of curvature A is a horizontal lift for the submersion 
T : M(k) —> N(k) of a curve with geodesic curvature — 2A in N(k), see [55, Lem. 3.1] and [44, 
Thm. 1.26]. In this paper we will only need the following facts about CC-geodesics in M(«) that 
can be proved without using their explicit expressions. 

Proposition 2.5. Let 7 be a complete CC-geodesic in M(re). Then we have 

(i) If n ^ 0 then 7 leaves any compact set of M(k) in finite time. 

(ii) If n > 0 then 7 cannot be contained inside an open hemisphere of § 3 . 

(iii) If k > 0 and 7 is closed, then the length of 7 is greater than or equal to 27r/ t/k. Moreover, 
equality holds if and only if 7 is a horizontal great circle of § 3 . 

Proof. For any k / 0 we can apply Lemma 2.4 (ii) to get that a CC-geodesic of curvature A and 
length L in M(k) is a reparameterization of a CC-geodesic of curvature A/-\Z[k[ and length L 
in M(k/|k|). So, we only have to prove the claim for k = —1,0,1. After a translation in the Lie 
group (M(k), *) we can suppose that 7(0) is the identity element o. By changing the orientation of 
7 we can admit that A ^ 0. 

We first treat the case n = —1,0. If A = 0 then 7 is a horizontal Riemannian geodesic, i.e., a 
straight line containing o inside the plane t = 0. If A ^ 0 then we know from [55, Lem. 3.1] that 
7 is a horizontal lift for the submersion T : M(k) —> N(k) of an arc c\ of a circle containing o. 
This arc meets <9 N(k) when k = — 1 and |A| < 1. For any z £ c\, z ^ o, let fL C N(«) be the 
circular sector bounded by the arc oz and the segment Hz. Then, the horizontal lift of oz through 
o is contained in 7, and we deduce by Lemma 2.3 that the f-coordinate of the point of 7 projecting 
onto z is proportional to the area enclosed by Q z . Statement (i) then follows since this area tends 
to 00 when the arc-length of oz goes to 00. 

Suppose now that k = 1. Again by [55, Lem. 3.1] we can find a circle c\ in N(l), bounding a 
disk of area Aq ^ 7r/2, and such that 7 is the horizontal lift of c\ through o. Let F be the oriented 
fiber of F containing o. From Lemma 2.3 the intersection 7 n T is a family of points at a constant 
distance 2Ao ^ n. Assume that 7 is contained in an open hemisphere E + associated to a linear 
hyperplane E of R 4 . It is easy to check that TO-E consists of two different points at distance 7r along 
r. Hence, it would be possible to find p £ yflT such that p £ E~ , a contradiction. This proves (ii). 
Finally, note that 7 is closed if and only if there are n,m £ N such that uAq = mir. From here we 
get n Js 2m since Ao ^ 7r/2. Moreover, we can take n in such a way that the lengths L and Lq of 7 
and c\ satisfy L = uLq. Hence the isoperimetric inequality in the surface N(l), which is a 2-sphere 
of curvature 4, implies 

L 2 = 71 2 Lq ^ n 2 (47tAo — 4Ag) = 47 t 2 to (n — m) ^ 47r 2 m 2 ^ 4-7T 2 , 
with equality if and only if m — 1, n = 2 and Aq = 7r/2, i.e., 7 is a horizontal great circle. □ 

It was proved in [54, Lem. 3.5] that the infinitesimal vector held associated to a C 2 family of 
CC-geodesics with the same curvature in the Heisenberg group M(0) satisfies a second order ODE 
analogous to the Jacobi equation in Riemannian geometry. The following result in [55, Lem. 3.3] 
establishes the same fact under weaker regularity conditions in arbitrary Sasakian 3-manifolds. 

Lemma 2.6. Let M be a Sasakian sub-Riemannian 3-manifold. We consider a C 1 curve a : I —> M 
defined on some open interval / C R, and a C 1 unit horizontal vector field U along a. For fixed 


PANSU SPHERES IN SUB-RIEMANNIAN 3-SPACE FORMS 


9 


A £ 1, suppose that we have a well-defined map F : I x I' M given by F(e,s) := 7 £ (s), where 
I' is an open interval containing the origin and 7 £ (s) is the CC-geodesic of curvature A with initial 
conditions 7 e ( 0 ) = a(e) and 7 e ( 0 ) = 17(e). Then, the vector field T4(s) := (dF/de)(e, s) satisfies 
the following properties: 

(i) V e is C°° along 7 e with [ 7 e , V e \ = 0, 

(ii) the function A (V e ,T^ + (14, 7 e ) is constant along r y e , 

(iii) V e satisfies the CC-Jacobi equation 

v" + R(%, V e )% + 2 A (j(V‘') - (V e , 7 e) T) = 0 , 

where R is the curvature tensor in (M, g) and the prime ' denotes the covariant derivative 
along the geodesic 7 e , 

(iv) the vertical component of V e satisfies 

(2.9) (V e ,T)' = 2 (V e ,J(%)), 

(2.10) (V e ,T)"' + 4 (A 2 + K) (V e , T)' = 0, 

where K is the Webster scalar curvature of M and ' is the derivative with respect to s. 

Remark 2.7. Following [34, Lem. 2.1] we can prove similar properties for the vector field 14 (.s) := 
(dF/de)(e,s) associated to a variation F(e,s) := 7 e (s), where 7 e (s) is the Riemannian geodesic in 
( M,g ) with 7 e ( 0 ) = a(e) and 7 e ( 0 ) = U(e). In this case, the second order equation in Lemma 2.6 
(iii) becomes the classical Jacobi equation 

(2.11) Vf' + R( ie,V e )% = 0. 


If 7 is a CC-geodesic in M then we define a CC-Jacobi field along 7 as a solution to the CC-Jacobi 
equation in Lemma 2.6 (iii). If the Webster scalar curvature of M is constant along 7 then, an easy 
integration from ( 2 . 10 ), gives us the following result. 


Lemma 2 . 8 . Let 7 (s) be a CC-geodesic of curvature X in a Sasakian sub-Riemannian 3-manifold 
M. Let V be the CC-Jacobi field associated to a variation ofj by CC-geodesics of curvature A as in 
Lemma 2.6. Suppose that the Webster scalar curvature of M is a constant n along 7 . If we denote 
t := 4 (A 2 + k) and v(s) := (V, T)(s), then we have 


(i) For t < 0, 


1 


v(s) = _ (a sinh(v / —rs) + b cosh(v / —r s)) + c, 


T 


where a := i/( 0 ), b := (1/y/^r) v"{0) and c := u( 0 ) + ( 1 /r) u"( 0 ). 

(ii) For t = 0, 

v (s) = as 2 + bs + c, 

where a = (1/2) v"(0), b = i/( 0 ) and c = u( 0 ). 

(iii) For r > 0, 

v(s) = —.= (a sin(v / rs) — b cos(v / rs)) + c, 

Vt 

where a = v'(ff), b = ( 1 /i/r) v"(0) and c = v( 0 ) + (1 /r) u"( 0 ). 


Remark 2.9. Previous work on the CC-geodesics equation, the second variation of length and CC- 
Jacobi fields in contact sub-Riemannian 3-manifolds is found in [33] and [56]. In [9], Chanillo and 
Yang studied the CC-Jacobi fields associated to C 2 families of CC-geodesics of the same curvature 
leaving from the same point. 
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2.4. Surfaces and singular set. 

Let M be a contact sub-Riemannian 3-manifold and E a C 1 surface immersed in M. Un¬ 
less explicitly stated we shall always assume that 9E = 0. The singular set of E is the set 
Eo := {p £ E; T p E = Tip}. Since LL is a completely nonintegrablc distribution, it follows by 
the Frobenius theorem that Eo is closed and has empty interior in E. Hence the regular set E — Eo 
of E is open and dense in E. By using the arguments in [19, Lem. 1], see also [2, Thm. 1.2] and [31, 
App. A], we can deduce that, for a C 2 surface E, the Hausdorff dimension of E 0 in ( M,g ) is less 
than or equal to 1. In particular, the Riemannian area of Eo vanishes. 

If E is orientable and N is a unit vector normal to E in (M, g), then we define the area of E by 

(2.12) A(E) := f \N h \dS, 

Jy, 

where dE is the Riemannian area element. If E bounds a set Q C M then the Riemannian diver¬ 
gence theorem implies that A(E) coincides with the sub-Riemannian perimeter a la De Giorgi of 
12, which can be introduced by following [7], see equation (6.1) for a precise definition. It is clear 
that So = {p £ E; Nh(p) = 0}, where Nh := N — (N, T)T. Hence, the differentiation of the 
area functional could be problematic when Eo 0. In the regular part E — Eo, we can define the 
horizontal Gauss map Vh and the characteristic vector field Z , by 

(2.13) v h := Z := J{v h ). 

As Z is horizontal and orthogonal to Vh : then Z is tangent to E. Hence Z v generates T p E D Tip for 
any E — Eo- The integral curves of Z in E — Eo are called ( oriented ) characteristic curves of E. 
They are both tangent to E and to H. If we define 

(2.14) S:=(N,T)vH-\N h \T, 

then {Z p , S p } is an orthonormal basis of T p E whenever p £ E — Eo- Moreover, for any p £ E — Eo 
we have the orthonormal basis of T p M given by { Z v , (u/j) p , T p }. From here we deduce the following 
identities on E — Eo 

(2.15) \N h \ 2 + (N,T) 2 = 1, {v h ) T = (N,T) S, T r = -\N h \S , 

where U T stands for the projection of a vector field U onto the tangent plane to E. 

If E is an oriented C 2 surface immersed in M then, for any p £ E — Eo and v £ T p M , the following 
equalities hold (see [34, Lem. 3.5] and [55, Lem. 4.2]): 

(2.16) v (|iV fc |) = (D V N, vh) + (N, T) (v , Z), 

(2.17) v((N, T» = (. D V N, T) + (N, J(«)>, 

(2.18) D v v h = llV^r 1 ((DyN, Z) - (N, T) (v, v h )) Z + (v, Z) T. 

We denote by B the shape operator of E in (M,g). It is defined, for any vector U tangent to E, 
by B(U) := —DuN. We will say that E is complete if it is complete with respect to the Riemannian 
metric induced by (M,g). 

2.5. Differentiation under the integral sign. 

In the next lemma we recall the classical Leibniz’s rule for differentiating the integral of a function 
which depends on a real parameter. This will be used in Sections 4 and 5 to compute the deriva¬ 
tives of the sub-Riemannian volume and area functionals. A proof can be obtained by following the 
arguments in [22, Thm. 2.27]. 

Lemma 2.10. Let (f2,/x) be a measure space and I C R an open interval. Consider a real valued 
function / :/ xll-> R such that the following conditions hold: 
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(i) for any s £ I, the function f(s , •) : S 2 —> R is integrable, 

(ii) the function f{-,x) : I —> R. is absolutely continuous for almost every x £ Q, 

(iii) there is an integrable function h : Q -7 R such that \(df /ds)(s,x)\ ^ h(x) for almost every 
s £ I and almost every x £ fi, 

(iv) there is so £ I such that (df /ds)(so, x) exists for almost every x £ fb 

Then, for the function A : I —> K defined by H(s) := f n f(s,x) dp(x), the derivative -A'(so) exists, 
and we have H'(so) = J a (df /ds)(so, x) dp{x). 

Example 2.11. Here we show that Lemma 2.10 need not hold if condition (ii) fails. Consider the 
function / : (—1,1) x (—1,1) —> R given by 


f{s,x) := 


— 1, — 1 < s < x, 
1, X < s < 1. 


Clearly H(s) = J_ x f(s, x) dx = 2s, so that T'(s) = 2 for any s £ (—1,1). Note that the conditions 
(i), (iii) and (iv) in Lemma 2.10 are satisfied. However, differentiation under the integral sign fails; 
otherwise we would obtain H'(s) = 0 for any s £ (— 1 , 1 ), a contradiction. 


3. The spherical surfaces S\(p) 

Pansu conjectured in [48, p. 172] that any isoperimetric surface in the first Heisenberg group M(0) 
is congruent to a sphere of revolution S\ obtained by vertical rotations of a CC-geodesic of curvature 
A ^ 0 connecting two points p and p\ of the vertical axis. Since vertical rotations are isometries 
of M(0), it follows from Lemma 2.4 (i) that S\ coincides with the union of all the CC-geodesics 
of curvature A joining p and p\. This approach was also employed in [35, Sect. 5.1] to produce 
spherical surfaces in the sub-Riemannian 3-sphere M(l). In this section we extend the construction 
inside any 3-dimensional space form M. For that, we must study first when the CC-geodesics of a 
fixed curvature A £ I leaving from a given point p £ M meet again at another point p\ £ M. 

Lemma 3.1. Let M be a 3-dimensional space form of Webster scalar curvature n. Fix a point 
p £ M and a number A £ K. Then, all the CC-geodesics of curvature A leaving from p meet again in 
M if and only if A 2 + n > 0. In such a case, the length of any CC-geodesic of curvature A between 
p and the first cut point p\ equals n/V A 2 + k. 

Proof. Fix a positive orthonormal basis {ei,e 2 } in the contact plane Up. For any 6 £ R, let 
70 : R — > M be the CC-geodesic of curvature A with 70 (0) = p and 70 ( 0 ) = U(9), where 
U{9) := (cos0)ei + ( 81110 ) 62 . The C°° map F : R 2 — > M defined by F(6,s) := 70 (s) is the 
flow of the CC-geodesics of curvature A leaving from p. We denote Vg(s) := (8F/d9)(9, s). This is a 
CC-Jacobi held along 7 g by Lemma 2.6 (iii). Let vg := (Ve,T}. By using Lemma 2.6 (ii), equation 
(2.9) and that 1^(0) = 0, we deduce 

(3.1) Vg = —(Avg) 70 + (v'g/2) J{ 70 ) + VgT. 

Evaluating at s = 0 we get Ug(0) = 0. We compute Vy (0) from (3.1), so that we obtain 

K (0 ) = ^J{U(9)). 

On the other hand, Vg = Dj e Vg = Dy g 70 since [ 70 , Vg] = 0 along 70 . Therefore 

70 (0) = - (sin 9) ei + (cos 9) e 2 = J (U (0)), 

0 
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so that Vg (0) = 2. Since M has constant Webster scalar curvature, we apply Lemma 2.8 to get 

f 2 (A^+«) (- 1 ~ cosh ( 2 + K ) s ))i if A 2 + /t<0, 

(3.2) vg(s) = ^ s 2 , if A 2 + k = 0, 

1 2 fA^y I 1 “ cos ( 2 a/A^+ks)), if A 2 + k > 0. 

It follows that vg does not depend on 9 , and that vg(s) ^ 0 for any s ^ 0 if A 2 + k ^ 0. For A 2 + 

k > 0 it is clear that the first so > 0 for which vg(so) = 0 is given by so := ir/y/X 2 + k. 

Let us prove the statement. If A 2 + k > 0 then (3.1) implies V^(sq) = 0 for any 0 £ R since 

vg(sg) = v' g (s o) = 0. This is equivalent to the existence of a first point p\ G M such that jg(so) = p\ 
for any 9 G K. Conversely, suppose that all the CC-geodesics jg meet again at p\ G M. This means 
that there is a C 1 function s(d) > 0 such that p\ = 7 e(s(0)) = F(6 , s{9)) for any 9 € R. Note that 
s(9) is the length of yg between p and p\. By taking derivatives with respect to 9 , we obtain 

O = y,(s(0))+s'(0)7s(s(0)), 

and so the vertical component vg of Vg vanishes at s(ff). From the explicit computation of vg in 

(3.2) , we conclude that A 2 + n > 0 and s{9) must be constant. Furthermore, the first cut point of 

the CC-geodesics ye is attained for s 0 = t/VA 2 + k. □ 

The previous lemma motivates the following definition. Let M be a 3-dimensional space form of 
Webster scalar curvature k. For any p £ M and any A ^ 0 such that A 2 + k > 0, we define the 
spherical surface S\(p) as the set 

(3.3) S\(p) := { 7 „(s); v G Hp, |n| = 1, 0 < s < tt/\ / X 2 + k}, 

where 7 „ : M —> M is the CC-geodesic of curvature A with 7 „( 0 ) = p and 7 „( 0 ) = v. The 
south and north poles of S\(p) are respectively given by p and p\ := 7 „( 7 t/v / A 2 + k) (which does 
not depend on v G T~L P with |t>| = 1). For further reference we define the equator of S\(p) by 
{7i)( 7r /(2 VX 2 + k)) ; v G hip, |i?| = 1} and the open hemispheres as the sets 

S\(p)~ := ( 7 „(s); v G Up, |i>| = 1, 0 < s < 7 r /(2 \/A 2 + k )}, 

S\{p) + ■= {y«(s) ; v G Tip, H = 1, 7 r /(2 \/A 2 + k) < s ^ t/v 7 A 2 + k}. 

As an immediate consequence of Lemma 2.4 we can deduce the behaviour of S\(p) with respect 
to local isometries and the homothetic deformation in ( 2 . 6 ). 

Lemma 3.2. Let M be a 3-dimensional space form of Webster scalar curvature n. For any p G M 
and any A ^ 0 such that A 2 + k > 0, we have 

(i) If 4>: M —► M' is a local isometry then < j)(S\(p )) = S' x {(f>{p)). 

(ii) For any e > 0, the set S\(p) associated to the sub-Riemannian metric gh coincides with the 
set S x/e (p) for the sub-Riemannian metric (gh)e := £ 2 9h- 

Our main goal in this section is to prove some properties of the sets S\ (p) regarding their sym¬ 
metry, embeddeness, topology and smoothness. They are gathered in the next result. 

Theorem 3.3. Let M be a 3-dimensional space form of Webster scalar curvature k. Fix a point 
p G M and a number A ^ 0 such that A 2 + k > 0. Then we have 

(i) The set S\(p) is invariant under any (f> G Iso(M) with 4>(jp) = p. Moreover, if M is homo¬ 
geneous, then any two sets S\(jp) and S\(q) are congruent. 

(ii) If M = M(k) then S\(p) is an embedded C 2 surface congruent to a sphere of revolution with 
respect to the vertical axis. 

(iii) The set S\(p) is a C 2 sphere immersed in M. This sphere is also C°° off of the singular 
set, which consists of the south and north poles. 
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(iv) The sphere S\(p) is always Alexandrov embedded. Moreover, there is Ao = Ao (p) ^ 0 such 
that Aq + k > 0 and S\(p) is embedded for any A > Ao- 

Remark 3.4. A compact surface E immersed in M is Alexandrov embedded if we can extend the 
immersion : E — > M to an immersion <po : T> —s > M of a compact 3-manifold T> with dT> = E. 


Proof of Theorem 3.3. Statement (i) is an immediate consequence of Lemma 3.2 (i). By using that 
M(k) is homogeneous and Lemma 3.2 (ii) it suffices to prove (ii) when n = —1, 0,1 and p equals the 
identity element o in (M(k),*). We shall denote the associated set by <5>£(o). 

The cases k = 0,1 are well known, see for example [54, Ex. 3.3] and [35, Prop. 5.1, Thm. 6.4]. 
For the remaining case k = — 1 we proceed as follows. Fix A > 1. Let 7 be a CC-geodesic in M(— 1) 
of curvature A and length n/x/X 1 — 1 leaving from o. By using [55, Lem. 3.1] or [44, Thm. 1.26] we 
deduce that the Euclidean projection of 7 onto the hyperbolic plane N(—1) is a clockwise parame¬ 
terized circle c\ of radius 1/(2A). Take a point q := ( z,t) G 7. Let f l~ C N(— 1) be the circular sector 
bounded by the arc oz of c\ and the segment 02 . By Lemma 2.3 we get that t/2 equals the area of 
in N(— 1). To compute this area as a function of |.z| we can use the Gauss-Bonnet formula. It follows 
that iS^(o) is the union of two radial graphs over the horizontal disk D\ := [z£l 2 ; |z| ^ 1/A}. 
In fact, a suitable vertical translation maps Sf(o) into {(z,t) G D\ x R.; t = ±f(\z\)}, where 
/ : [0,1/A] —> R is the function given by 


/(*) := o 1- 


A 


A 




aresm 


\/\ 2 — 1 x 
\/l — x 2 


— arctan 


A: 


x/1 — X 2 x 2 


From here it is easy to check that S£ (o) is a C 2 embedded sphere of revolution with singular points 
only at (0, ±/(0)). Moreover, <S^(o) is C°° off of these points, where it fails to be C 3 . This is also 
known in M(0), see the remark after Thm. 8.16 in [ 8 ], and in M(l) for A ^ 0, see [35, Rem. 5.2]. 


Now we can prove (iii) and (iv). We apply Proposition 2.1 to the 3-dimensional space form M, 
so that we find a surjective local isometry II : M(re) —>• M. Given a point po G M(/c) with II(po) = P , 
we have by Lemma 3.2 (i) that S\(p) = Il(<S£(po))- This allows to deduce (iii) from (ii). Moreover, 
S\(p) is Alexandrov embedded since the sphere S£(po) bounds a topological 3-ball in M(re). Finally, 
the fact that II : M(k) — > M is a local isometry implies the existence of an open ball B := B K (po, Rq) 
for the Carnot-Caratheodory distance in M(k) such that II : 23 —>- 11(1?) is an isometry. We denote 
by Ao the unique non-negative number with Ag + k > 0 and n / yj Aq + k = Rq. By definition of 
CC-distance it is clear that Sf(po) C B for any A > Ao. In particular, the sphere S\(p) = n(<S£(po)) 
is embedded for any A > Ag. This completes the proof. □ 


Remarks 3.5. 1. In a homogeneous Sasakian 3-manifold M the number Ao = Ao(p) in statement 
(iv) does not depend on p. So, the spheres S\(p) are embedded for any A > Ao and any p G M. 

2. If M has non-trivial topology, then the spheres S\(p) need not be embedded for all values of 
A. To illustrate this we consider a Sasakian cylinder M := M(k)/G as the ones in Example 2.2. 
Take a sphere of revolution Sf(o) in M(/c) with poles identified under the action of G. If we denote 
p := II(o) then S\(p) = 11 ( 5 ^ (o)), so that S\(p) is immersed in M but not embedded. 

3. The arguments in Theorem 3.3 allow also the construction of immersed planes inside any 
3-dimensional space form M. Given p G M and A > 0 such that A 2 + k ^ 0, we define the set 

(3.4) C x {p) ■■= { 7 u( s ) Up, |u| = 1 , s ^ 0 }, 

where j v : R —> M is the CC-geodesic of curvature A with 7 „( 0 ) = p and 7 „( 0 ) = v. By equation 
(3.2) we know that the vertical component of the CC-Jacobi field given in (3.1) never vanishes for 
s > 0. Hence, the set C\{jp) is a complete C 2 plane immersed in M whose singular set is {p}. 

4. It is natural to ask if our construction of S\(p) would provide spherical surfaces inside an 
arbitrary complete, contact, sub-Riemannian 3-manifold M of constant Webster scalar curvature. 
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Observe that, if M is not Sasakian, then the characteristic curves of a CMC surface are not neces¬ 
sarily CC-geodesics, but they still satisfy a second order ODE similar to (2.8), see [24, Prop. 4.1]. 
Hence we may consider the set S\(p) of all the solutions to such ODE leaving from p and study, as 
in Lemma 3.1, when these solutions meet again. However, the analysis of the associated CC-Jacobi 
fields is much more involved since their vertical component does not only depend on the Webster 
scalar curvature but also on the Webster-Tanaka torsion [24, Prop. 4.3]. 


We conclude this section by introducing polar coordinates inside the spherical surfaces S\ (p) and 
giving some computations in such coordinates that will be helpful in the sequel. 

Lemma 3.6. Let M be a 3- dimensional space form of Webster scalar curvature k. Given a point 
p £ M, a positive orthonormal basis {ei, 62 } in Tip, and a number A Js 0 with A 2 + k > 0, we define 
the map F : R 2 — > M by F(6, s ) := 70 (s), where 70 is the CC-geodesic of curvature A with 70 ( 0 ) = p 
and 70 ( 0 ) = (cos (9) e\ + (sin0) e 2 - Denote t := V\ 2 + k. The following properties are satisfied: 


(i) the restriction of F to C\ := [0,2-7r] x (0,7t/t) is a C°° immersion whose range is S\(p) 
minus the poles, 

(ii) there is a Riemannian unit normal N along S\ (p) such that any CC-geodesic 70 (s) with 
s £ (0, 7 r/r) is a characteristic curve of S\{p), 

(iii) with respect to the coordinates ( 9,s ) £ C\, these equalities hold 


dS\{p) = /1 _|_ (^_2 — 1 ) cos 2 ( rs ) dOds, 


\N h \ = 


sin(rs) 


y/l + (t 2 — 1) COS 2 (ts) 

(B(Z), Z) = 2A \N h \, (B(Z ), S) = (1 - r 2 ) \N h \ 2 , 


(N,T) 
< B(S),S ) 


T cos(rs) 

\/l + (t 2 — 1) cos 2 (ts) ’ 
A t 2 \Nh\ 

1 + (r 2 — 1) cos 2 (rs) ’ 


where dS\(p) is the Riemannian area element, {Z, S'} is the orthonormal basis in (2.13) and 
(2.14), and B is the Riemannian shape operator with respect to N. 


Proof. The map F : R 2 — > M is C°° since, for fixed A, the solutions of (2.8) depends differentiably 
( C°° ) on the initial data. Note that ( 1 dF/ds){0, s) = 70 (s), which is a horizontal vector. On the 
other hand, Vg(s) := ( dF/dd)(6,s ) is a CC-Jacobi field along 70 by Lemma 2.6 (iii). From the 
computations in the proof of Lemma 3.1 we know that 

Vg = -(Aw) 70 + [v'/ 2) J( 70) + vT, 


where the prime ' denotes the derivative with respect to s, and v is the function given by 


(3.5) 


;(s) := 


1 — cos( 2 rs) sin 2 (ts) 


2 T 2 T 2 

The fact that v(s) > 0 for s £ (0,7t/t) implies that the differential of F has 
(0, s) £ C\, and so F : C\ —> M is a C°° immersion. By the definition of S\{p) 
that F(C\) = S\(p) — {p,p\}. This proves (i). Now, we define 


rank two for any 
in (3.3) it is clear 


(3.6) 


-wJ(70) + K/2)T 
\Jv 2 + (w'/2) 2 


which provides a Riemannian unit normal along S\(p) — {p,p\}. By using (2.13) we get Vh = —*^( 70 ) 
and Z = 70 . This means that any 70 (s) with s £ (0,7t/t) is a characteristic curve of S\(p), proving 
statement (ii). On the other hand, it is clear that 


dS x (p) = y/\V e \ 2 - \ 2 v 2 dOds = -y/w 2 + (w'/2) 2 ddds 

= ./ sin 2 (ts) + t 2 cos 2 (ts) d6 ds = 1 + (t 2 — 1 ) cos 2 (ts) dO ds. 
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Moreover, the expressions of \Nh\ and (IV, T) in statement (iii) can be obtained from (3.6) and (3.5). 
Note that 0 = Z((N, Z)) = (J D Z N, Z) + (iV, DzZ ). From equation (2.8), we get 

(B(Z), Z) = (TV, D z Z ) = —2A (TV, J(Z)) = 2A (TV, v h ) = 2A |TV h |. 

By taking into account (2.17), we obtain 

Z((N,T)) = < D Z N,T) + (N,D Z T) = \N h \ (B(Z),S) - |TV fc |, 

S((N, T)) = ( D S N , T) + (TV, £ S T> = |JV fc | (B(S), S ), 
since the tangent projection of T is — \Nh\ S, D Z T = —Uh and DsT = (TV, T) Z. It follows that 


' ( h } \N h \ + |JV fc | + 1 + (r 2 - l)cos 2 (rs) 


(B(S), S ) = S( ^’P ) = -A Z((N, T» = -A (TV, T)' = 


= (l-r 2 )|iV,| 2 , 
At 2 sin(rs) 


1^1 (l + ( r 2 _i) cos 2 (ts )) 3/2 ’ 

where in the second equation we have used that (TV, T) does not depend on 9 , and that 


(3.7) 

This completes the proof. 


S = 


- 1 


VW+(« 


: - A |JV„| Z. 


□ 


An easy consequence of Lemma 3.6 is the following fact, that will play a key role in this work. 

Corollary 3.7. Let S\(p) be a spherical surface inside a 3 -dimensional space form. Then, for any 
Riemannian unit normal vector N along S\(p), the function |7V/j| _1 is integrable along S\(p) with 
respect to the Riemannian area element dS\ (p). 


4. Characterization results and classification of stationary surfaces 

Let M be a 3-dimensional space form. Our main goal in this section is to establish characterization 
results for the spherical surfaces S\(p) introduced in (3.3) as critical points of the sub-Riemannian 
area under a volume constraint. More precisely, we will be able to prove sub-Riemannian coun¬ 
terparts to the classical theorems of Hopf and Alexandrov for constant mean curvature surfaces 
in Riemannian 3-space forms. We will also provide the description of complete volume-preserving 
area-stationary surfaces in M with non-empty singular set, thus extending previous results in [54] 
and [35] . We first introduce basic definitions and review some properties of stationary surfaces that 
will be key ingredients throughout this section. 

4.1. Volume-preserving area-stationary surfaces. 

Let M be a Sasakian sub-Riemannian 3-manifold and tpo : E — > M an oriented C 2 surface im¬ 
mersed in M. By a variation of £ we mean a map ip : I x £ M (which is assumed to be C 2 
unless otherwise stated), where / C R is an open interval containing the origin, and ip satisfies: 

(i) tp{0,p) = <p 0 {p) for any pe£, 

(ii) the map ip s : £ — > M given by tp s {p) '■= (p(s,p) is an immersion for any s € /, 

(iii) there is a compact set C C £ such that p s (p) = tpo(p) for any s £ I and any p £ £ — C. 

The velocity vector field associated to ip is the C 1 vector field along £ defined by U p := (d<p/ds)(0,p). 
The area functional is given by A(s) := A(£ s ) = A(y> s (£)). We define the volume functional V(s) 
as in [3, Sect. 2], by taking the signed volume enclosed between £ and £ s . More precisely 

(4.1) V(s) := [ ip*(dM), 

J[0,s]xC 
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where cLM is the Riemannian volume element in (M,g). The variation is volume preserving if 
V (s) is constant for any s small enough. We say that E is volume-preserving area-stationary or 
area-stationary under a volume constraint if A' (f)) = 0 for any volume-preserving variation of E. 

Let us recall the computation of A'(0) and V'(0). From (2.12) and the fact that Eo has vanishing 
Riemannian area, we have 

(4.2) -A(s) = [ |-/W| p (s) |JacvJ s |pdE. 

J E-Eo 

Here |iV/i|p(s) := |iV) l |(<p s (p)), where N is a C 1 vector field along the variation whose restriction to 
any E s provides the Riemannian unit normal N s compatible with the orientations of E and M. If 
p £ E and {ei,e 2 } is any orthonormal basis in T p Y, then |Jac</j s | p denotes the squared root of the 
determinant of the positive definite matrix G(s) := ((e»(y? s ), e J -(^ s ))) i . with i,j = 1,2. Note that 
\N h \ p (s) is a positive Lipschitz function, for s small enough. Moreover, we have 

l-^lp( s ) = £f(l-Nh|) = ( DuN h ,u h ), |Jac(p s |p(s) = i i/det G(s) trace(G'(s) G(s) _1 ), 

where Uh is the horizontal Gauss map defined in (2.13). The previous expressions are uniformly 
bounded as functions of s and p near s = 0. Hence we can use Lemma 2.10 to differentiate under 
the integral sign in (4.2). By reproducing the arguments in [54, Lem. 4.3], we get 

(4.3) A'(0) = —2 f H(U,N)dE- f div E «£/, N) K) T ) dE, 

—So J S—So 

provided the function 

(4.4) - 2 H := div s Vh 

is locally integrable with respect to dE. Here div E Vh and (vh) T are the Riemannian divergence 
relative to E and the tangent projection of Vh, respectively. The function H in (4.4) is the ( sub- 
Riemannian) mean curvature of E, as defined in [53], [34] and [55]. On the other hand, it is known 
[3, Lem. 2.1] that V'(0) is given by 

(4.5) V"(0) = f (U,N)dS. 

Jy. 

We say that E is a constant mean curvature surface (CMC surface) if the function H in (4.4) is 
constant on E — Eo- As an easy consequence of (4.3) and (4.5) it follows that E is volume-preserving 
area-stationary if and only if there is a constant Asl such that (A + 2AH)'(0) = 0 for any variation 
of E. In particular E has constant mean curvature A, see [55, Prop. 4.3] for details. 

The regular part E — Eo of a CMC surface is ruled by CC-geodesics of M, see for example [54], 
[13] and [30]. In precise terms we get this result, which can be proved as in [54, Thm. 4.8]. 

Proposition 4.1. Let E be an oriented C 2 surface immersed inside a Sasakian sub-Riemannian 
3-manifold. Then, in the regular part E —Eo, we have equality DzZ = (2 H)vh- Therefore, ifT, has 
constant mean curvature H, then any characteristic curve of E is a CC-geodesic of curvature H. 

The singular set E 0 of a CMC surface E is well understood by the results in [13] for arbitrary 
pseudo-Hermitian 3-manifolds, see also [24, Sect. 5]. By using the ruling property in Proposition 4.1, 
the description of E 0 given in [13] can be stated as follows. 

Theorem 4.2 ([13, Thm. B]). Let E be an oriented C 2 surface of constant mean curvature H im¬ 
mersed in a Sasakian sub-Riemannian 3-manifold M. Then, the singular set Eo consists of isolated 
points and C l curves with non-vanishing tangent vector (singular curves). Moreover, we have: 
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(i) ([13, Thm. 3.10 and Lem. 3.8]) if p £ Eo is isolated, then there exist r > 0 and A £ M with 
|A| = \H\ such that the set D r (p) := { 7 „(s); v £ TL P , |u| = 1, s £ [0,r)} is an open neigh¬ 
borhood of p in E. Here j v denotes the CC-geodesic of curvature A in M with 7„(0) = p 
and 7„(0) = v. Moreover, the index of the vector field |7V^| Z at p equals +1. 

(ii) ([13, Prop. 3.5 and Cor. 3.6]) if p is contained in a C 1 curve P C Eo then there is a 
neighborhood D of p in E such that D is the union of two disjoint connected open sets 
D + ,D~ C E —Eo- Moreover, for any q GTflD there are exactly two CC-geodesics 71 C D + 
and 72 C D~ of curvature A leaving from q and meeting transversally P at q with opposite 
initial velocities. The curvature A does not depend on q and satisfies |A| = \H\. 

Remarks 4.3. 1. The relation between A and H depends on N p . If N p = T p then A = H, and the 
CC-geodesics in the statement are characteristic curves of E. If N p = — T p then A = —H. 

2. The regular and singular sets of a C 1 surface in the Heisenberg group M(0) satisfying the 
constant mean curvature equation in a weak sense have been studied in [16] and [14], see also [27]. 

The previous result together with the Hopf index theorem for line fields implies the following 
topological restriction proved in [13, Thm. E], see also [24, Thm. 5.7]. 

Theorem 4.4. Let E be a compact, connected, oriented C 2 surface immersed in a Sasakian sub- 
Riemannian 3-manifold. IfT, is a CMC surface, then the genus ofT, equals 0 or 1. Moreover, E is 
homeomorphic to a sphere if and only if E contains at least one isolated singular point. 

In sub-Riemannian geometry a CMC surface need not be volume-preserving area-stationary. In¬ 
deed, we have the following characterization result, see [54, Thm. 4.17 and Prop. 4.20] for a proof in 
M(0) and [15], [12], [31], [24] for generalizations with low regularity assumptions and other settings. 

Theorem 4.5. Let E be an oriented C 2 surface immersed in a Sasakian sub-Riemannian 3-manifold. 
Then E is area-stationary under a volume constraint if and only if E has constant mean curvature 
and the characteristic curves meet orthogonally the singular curves when they exist. Moreover, in 
such a case, any singular curve in E is of class C 2 . 

Example 4.6. Let M be a 3-dimensional space form. Consider a spherical surface S\(p) as defined 
in (3.3). By using Lemma 3.6 (ii), equation (2.8) and Proposition 4.1 we deduce that S\(p) has 
constant mean curvature H = A. Moreover, as S\ (p) has two singular points we conclude from 
Theorem 4.5 that S\(p) is area-stationary under a volume constraint. The same argument shows 
that any plane C\(p) as in (3.4) is volume-preserving area-stationary with H = 0. 

4.2. Characterization results. 

Now we can prove a first characterization result for the spherical surfaces S\(p). This may be 
seen as a sub-Riemannian counterpart to Hopf uniqueness theorem, see [32, Chap. 6 ] and [17], which 
states that any constant mean curvature sphere immersed inside a complete Riemannian 3-manifold 
of constant sectional curvature must be round. 

Theorem 4.7 (A Hopf theorem for sub-Riemannian 3-space forms). If E is a topological sphere of 
class C 2 and with constant mean curvature immersed inside a 3-dimensional space form M, then E 
coincides with one of the surfaces S\(p) defined in (3.3). 

Proof. First we apply Theorem 4.4 to deduce that E contains an isolated singular point p. We 
choose the unit normal N to E such that N p = T p . Let n be the Webster scalar curvature of M, 
and A the mean curvature of E with respect to N. By Theorem 4.2 (i) and Remarks 4.3 we can 
parameterize a small open neighborhood of p in E by means of the flow of all the CC-geodesics 7 v of 
curvature A in M leaving from p. By completeness any 7 „ can be extended as a characteristic curve 
of E until it meets a singular point of E. If A 2 + k ^ 0 then equation (3.2) would imply that the 
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vertical component of the associated CC-Jacobi fields never vanishes along E — {p}. Thus, the point 
p would be the unique singular point of E, a contradiction with the Hopf index theorem for line 
fields since E is topologically a sphere. Hence A 2 + k > 0. If A ^ 0, then it is clear that E = S\(p). 
In the case A < 0 we obtain E = S-\(p \), where p\ is the first cut point of the CC-geodesics j v . □ 

Remark 4.8. Recently Cheng, Chiu, Hwang and Yang have proved a Hopf theorem in the Heisen¬ 
berg group H" with n ^ 2 under the additional assumption that E is umbilic [10, Cor. A']. 

The ideas in the previous proof can be used to obtain the classification of CMC surfaces with 
isolated singular points in 3-dimensional space forms. 

Theorem 4.9. Let E be a complete, connected and orientable C 2 surface immersed inside a 3- 
dimensional space form M. If E has constant mean curvature and at least one isolated singular 
point, then E is either a spherical surface S\(p) as in (3.3), or a plane C\(p) as in (3.4). 

In [54, Thm. 6.8] and [35, Thm. 5.9] it was proved, by means of an explicit computation of 
the CC-geodesics, that any singular curve of a C 2 volume-preserving area-stationary surface in the 
Heisenberg group M(0) or in the unit sphere M(l) is a piece of a CC-geodesic. Here we are able 
to generalize this fact to arbitrary 3-dimensional space forms by using the explicit expression of the 
vertical component of a CC-Jacobi field given in Lemma 2.8. 

Theorem 4.10. Let E be a complete and orientable C 2 surface immersed inside a 3 -dimensional 
space form M. IfT, is volume-preserving area-stationary then any singular curve ofT, is a complete 
CC-geodesic in M. 

Proof. Let T : / —> E be a connected and relatively compact open piece of a singular curve in E. By 
Theorems 4.5 and 4.2 we can suppose that T is a C 2 horizontal curve parameterized by arc-length. 
This implies that T' = hJ(T ), where h : I —> R is the continuous function defined by h := (T',J{T)'). 
To prove the statement it suffices by (2.8) to show that h is a constant function. 

We take a unit normal N to E for which N = T along T. We know that E has constant mean 
curvature A with respect to N. Let F : I x [0, +00) ->Mbe the C 1 map defined by F(e, s) := 7 e (s), 
where 7 e : [0, +00) —> M is the CC-geodesic of curvature A with 7 e (0) = T(e) and 7 e (0) = J(T(e)). 
By using Theorem 4.2 (ii), Remarks 4.3, Theorem 4.5 and the completeness of E, we deduce that 
any 7 e (s) with s > 0 is a characteristic curve of E until it meets a singular point. We denote 
IT(.s) := (dF/de)(e,s). From Lemma 2.6 (iii) this is a CC-Jacobi field along 7 £ . Let v e := ( V e ,T ). 
By Lemma 2.6 (ii), equation (2.9) and equality 14(0) = T(e), we get 

(4.6) 14 = ~{\v e ) 7 e + (u'/2) J(%) + v e T. 

Claim (*). Suppose that there is an open interval I'CI and a positive function s(e) defined on 
/' such that u e (s(e)) = 0 and I4(s(£)) 7^ 0 for any e £ I’. Then s(e) is constant on I’. 

To prove the claim we proceed as follows. Consider the curve T 0 (e) := F(e, s(£)), e S Two tan¬ 
gent vectors to E along T 0 are given by ( dF/de)(e , s(e)) = V e (s(e)) and ( dF/ds)(e , s(e)) = 7 e (s(e)). 
These vectors are horizontal and orthogonal by (4.6). As a consequence, To(£) is a singular curve 
of E. Up to a reparameterization we can suppose by Theorem 4.5 that To is a C 2 curve. Since the 
differential of F : /' x [0, +00) E has rank two for any pair (e, s(e)) we deduce from the inverse 
function theorem that s(e) is a C 1 function. Moreover, by the definition of To(e) we have 

f 0 (e) = I4(s(e)) + s'(e) j e (s(e)). 

As E is volume-preserving area-stationary, we conclude from the orthogonality condition in Theo¬ 
rem 4.5 that To(£) J_7 e (s(e)) and so, s'(e) = 0 for any e £ This proves the claim. 
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Now we compute the initial conditions of v E . Clearly w £ (0) = 0 since 14(0) = f(e). Evaluating 
(4.6) at s = 0 we see that u'(0) = —2. From (4.6) we also obtain 

K(0) = (2A) 7 e (0) + (v'J (0)/2) J(%(0)) - J(%)'(0) -2 T 7t(0) . 

Note that J (%)' = (2A)7 e — T by (2.2) and (2.8). Therefore 


F £ '(0) = -«(0)/2)f( e )-T r(£) . 

On the other hand, Lemma 2.6 (i) implies [%, V E ] = 0, so that Vf = D^ e V s = P>v e f E along y E . Hence 
K(0) = D Ve J{t) = J(t'(e))-T m =-(h(e)t( E ) + T m ), 


where we have used (2.2) and f' = h J(t). The two previous equalities for V4(0) yield u"(0) = 2 h(e). 


Let k be the Webster scalar curvature of M and denote t := 4 (A 2 + k). From Lemma 2.8 we can 
obtain the explicit expression of v e depending on the sign of r. Suppose that r < 0. Then, we have 


,(*) = 


h{e) 


(cosh(V^-^s) — l) — sinh(\/—r s) 


It follows that v e {s ) = 0 for some s > 0 if and only if </>(s) = h(e)/y/—T, where </> : (0,+oo) —> K 
is defined by </>(s) := sinh(^/—r s)/(cosh (^/—t s) — 1). Note that </> is decreasing on (0,+oo) and 
<^(0,+oo) = (1,+oc). Hence, there is a positive solution of equation v E (s) = 0 if and only if 
h(e) > y/—T. If h = \J—t on I then the statement is proved. Otherwise, we can find £o € / such that 
h(e o) ^ y/—T. After reversing the orientation of T if necessary we can suppose that h(e o) > \[—t. 
Consider the set A := {e S I ; h(e) = h(e o)}, which is clearly closed in I. Let us see that A is also 
open in I, then proving that h is constant. Take £\ € A. By continuity there is an open interval 
/'Cl with £\ € /' and h > \J~r on I'. Thus, for any e £ I', there is a unique s(e) > 0 for which 
v e (s(e)) = 0. From the definition of s(e) we get v' E (s(e)) = 2 and so, I4(s(e)) = J('y E (s(e))) ^ 0 by 
(4.6). By applying the Claim (*) we conclude that s(e) is constant on I'. Hence h(e) is constant 
on /' and /' C A , as desired. This completes the proof when r < 0. If r > 0 or r = 0, then the 
function v e (s) is respectively given by h(e) s 2 — 2s or 2r -1 / 2 [4(e) r -1 / 2 (1 — cos( v 4 r s)) — sin(y / rs)], 
and we can reason as in the case r < 0 to finish the proof. □ 


Now we can establish a sub-Riemannian counterpart to a classical theorem of Alexandrov [1], 
which states that any compact and embedded constant mean curvature surface inside the Euclidean 
space, the hyperbolic space or an open hemisphere of § 3 must be a round sphere. The result was 
already proved for the Heisenberg group M(0) in [54, Thm. 6.10]. 

Theorem 4.11 (An Alexandrov theorem in M(k)). Let E be a compact, connected, orientable C 2 
surface immersed in some model space M(k) with k ^ 0 or inside an open hemisphere of M(k) with 
n > 0. //E is volume-preserving area-stationary then E coincides with one of the spherical surfaces 
S\(p) defined in (3.3). 

Proof. Suppose that E has empty singular set or that it contains a singular curve. In both cases, 
Proposition 4.1 or Theorem 4.10 would imply that there is a complete CC-geodesic of M(k) included 
in E. By Proposition 2.5 we get a contradiction with the compactness of E or with the fact that E is 
contained inside an open hemisphere. Hence, we conclude that E must contain an isolated singular 
point. In particular, E is topologically a sphere by Theorem 4.4. Moreover, E has constant mean 
curvature since it is volume-preserving area-stationary. We finish by using Theorem 4.7. □ 

Remark 4.12. In Theorem 4.11 the hypothesis that E is contained inside an open hemisphere of 
M(k) when n > 0 cannot be removed. In fact, any Clifford torus given by the product of two circles 
with suitable radii is a volume-preserving area-stationary surface in M(re), see [35, Ex. 4.6]. Other 
examples of embedded CMC tori with rotational symmetry in M(l) were found in [35, Rem. 6.6]. 
Of course, none of these tori is inside an open hemisphere of M(k). Note that Theorem 4.11 also 
fails in arbitrary 3-dimensional space forms of Webster scalar curvature n ^ 0. Let M \= M {n)/G 
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be the Sasakian cylinder defined in Example 2.2. It is known, see [55, Sect. 4.2], that any Euclidean 
vertical cylinder E about the vertical axis is a CMC surface in M(re) with empty singular set. In M 
the quotient surface E /G provides a volume-preserving area-stationary embedded torus. 

The results in Theorems 4.10, 4.5 and 4.2 also allow to describe and classify the complete volume¬ 
preserving area-stationary surfaces with singular curves in any 3-dimensional space form. By using 
the same ideas as in [54, Thm. 6.11] and [35, Thm. 5.9] we can prove the following result. 

Theorem 4.13. Let E be a complete, connected and orientable C 2 surface immersed in a 3-dimen¬ 
sional space form. If E is volume-preserving area-stationary and contains at least one singular curve, 
then E coincides with one of the surfaces C^pfT) defined in Example below. 

Example 4.14. We follow the construction in [54, Ex. 6.7] and [35, Ex. 5.8]. Let M be a 3- 
dimensional space form and T a complete CC-geodesic of curvature p. At any point of T there are 
two unit horizontal vectors orthogonal to T. Then, for any A £ R, we can consider two families 
of CC-geodesic rays of curvature A £ R leaving orthogonally from T. We extend these geodesics 
until they meet a new singular curve “parallel” to T. Now, we repeat the construction with the new 
singular curves and so on. The resulting set C M ,A(r) is a volume-preserving area-stationary surface 
of constant mean curvature A whose singular set consists of CC-geodesics of curvature p. 

Note that Theorems 4.9 and 4.13 provide the classification of complete volume-preserving area¬ 
stationary C 2 surfaces with non-empty singular set in any 3-dimensional space form. This generalizes 
previous results in [54, Sect. 6] for M(0) and in [35, Sect. 5] for M(l). Some characterizations of 
CMC surfaces with empty singular set in M(re) can be found in [11], [4], [53], [54], [35] and [55]. 


5. Stability properties of S\(p) 

In this section we show that the spherical surfaces S\(p) introduced in (3.3) are second order 
minima of the sub-Riemannian area under certain deformations. We first state a second variation 
formula for S\(p) and discuss the technical problems arising in the computation. 

Let E be an oriented C 2 surface immersed inside a Sasakian sub-Riemannian 3-manifold M. For 
any variation p : I x E — > M we know from (4.2) that the associated area functional is given by 

A(s) = [ f(s,p)d,H, 

Je-e o 

where f(s,p) := |lW| p (s) |Jac</j s | p . In Section 4.1 we used differentiation under the integral sign, 
as stated in Lemma 2.10, to get A'(0) = f s _ s (df/ds)(0,p)dT,. This was possible since, for fixed 
p £ E — So, the map s H > f(s,p) is absolutely continuous, the derivative (df/ds)(0,p) exists, and 
( df/ds)(s,p ) is uniformly bounded as a function of s and p near s = 0. If we want to compute A w (0) 
from Lemma 2.10 we need some conditions that the derivatives (df /ds)(s,p) and (d 2 f / ds 2 )(s,p) 
must satisfy. This leads us to the notion of admissible variation that we now introduce. 

Definition 5.1. Let p : I x T, — I M be a C 1 variation of E. Recall that there is a compact set 
CCS such that p s (p) = 7>o(p) for any s £ I and any p £ E — C. We say that p is admissible if the 
function / : I x E —► R defined by f(s,p) := |lVh|p(s) \I&cp s \ p satisfies these properties: 

(i) (i df/ds)(s,p) exists for any s £ I and any p£E - Eo, 

(ii) the map s i—>■ (df/ds)(s,p) is absolutely continuous for almost every p £ E — Eo, 

(iii) there is a function h : E —> R, which is integrable on C with respect to dE, and such that 
|(<9 2 //9s 2 )(s,p)| ^ h(p) for almost every p € E - E 0 and almost every s £ I, 

(iv) the derivative (d 2 f /ds 2 )(0,p) exists for almost every p £ E — E 0 . 
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Clearly we can replace E — Eo with C C (E — Eo) i n the previous definition. Note that condition 
(i) holds provided the variation preserves the regular set of E, i.e., <p s (£ — £o) C £ s — (£ s ) 0 for 
any s £ I. In particular it is easy to check that, if E — Eo is a C 3 surface and tp is a C 3 variation 
compactly supported on E — Eo (which implies that tp leaves Eo invariant) then there is an interval 
/' CC I such that the restriction of tp to /' x E is admissible. In Appendix B we will show several 
examples of admissible variations for which Eo is not fixed. 

For a spherical surface S\(p) the family of admissible variations moving the poles is very large. 
In order to analyze if S\(p) is a second order minimum of the area we must obtain an explicit 
expression of A"{ 0) for all these variations. In [55, Thm. 5.2] this was done for variations compactly 
supported off of the poles and having the form p s {p) '■= exp p (u>(s,p) N p ). In the next result we 
extend the second variation formula to the case of arbitrary admissible variations. 

Theorem 5.2 (Second variation formula for S\(p)). Let M be a 3-dimensional space form of Web¬ 
ster scalar curvature k. Consider a spherical surface S\(p) for some p G M and A ^ 0 with 
A 2 + k > 0. Let p : I x S\(p) -A M be an admissible variation of class C 3 off of the poles. Denote 
u := ([/, N), where U is the velocity vector field and N is the unit normal in Lemma 3.6 (ii). Then, 
the functional A + 2XV is twice differentiable at s = 0, and we have 

(5.1) (A + 2AF)"(0) = T(u,u), 
where I is the quadratic form on C 1 (S\(p)) defined by 

(5.2) l(y,w) := [ \N h \~ 3 iz(v) Z(w) - (l A (r 2 - T)\N h \ 2 ) 2 vw\ dS\(p). 

Js x ( P ) 1 J 

In the previous formula Z is the characteristic field introduced in (2.13) and r := y/X 2 + k. 

That the expression in (5.1) coincides with the one obtained in [55, Thm. 5.2] is a consequence 
of Lemma 3.6 (iii), from which it is easy to check that 


(5.3) 


I B(Z) + S| 2 + 4 (« - 1) \N h \ 2 = (1 + (r 2 - 1) \N h \ 2 y 


along S\(p). Indeed, formula (5.1) also holds for some particular admissible variations with regular¬ 
ity less than C 3 off of the poles, see Remarks A.3 for precise references. Note that the integrability 
of lA^I -1 stated in Corollary 3.7 implies that l(v,w) is finite for any v,w € C 1 (S\(p)). We refer to 
I as the index form associated to S\(p) by analogy with the Riemannian situation studied in [3]. 

The proof of Theorem 5.2 is deduced from the more general second variation formula in Theo¬ 
rem A.l for admissible deformations possibly moving the singular set. After differentiation under 
the integral sign and a long computation, we get that (A + 2AF)"(0) equals I(u,u) adding the 
divergence of some tangent vector fields defined off of the poles. Then we use a generalized diver¬ 
gence theorem together with the integrability of |-/W| -1 to show that these extra terms vanish, thus 
proving (5.1). All the details are found in Appendix A. 

The next ingredient that we need is the expression of the index form with respect to the polar 
coordinates introduced in Lemma 3.6. 

Lemma 5.3. Let M be a 3-dimensional space form of Webster scalar curvature n. Consider a 
spherical surface <S\(p) for some p € M and A ^ 0 with A 2 + k > 0. For any function u £ C 1 (S\(p)) 
we denote u := uo F, where F(0 , s ) := 7 e(s) is the flow of CC-geodesics of curvature A leaving from 
p. Then, the index form of S\(p) defined in (5.2) satisfies 


1 


(£) 

where r := y/\ 2 + k and f(9,x) := \J\ + (r 2 — 1) cos 2 (a;) u(9,x/t). 


T{u,u) = - 

T J [0,27r] X [0,7r] 
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Proof. Let N be the unit normal along S\(p) for which any CC-geodesic 70 (s) with s £ (0,7t/t) 
is a characteristic curve of S\(p). From the definition of T(u,u) and the expressions for \Nh\ and 
dS\(p) in Lemma 3.6 (iii), we obtain 


(5.4) I(u,u) = 


' [0,27t] x [0,7r/r] 


1 + (t 2 — 1) cos 2 (ts) /i du\ 2 


V ds 


1 + (t 2 — 1 ) cos 2 (rs) 


u 2 


dO ds. 


We define the C 1 function u>(0, s ) := y/l + (r 2 — 1) cos 2 (ts) u(0, s). An easy computation gives us 
(5.5) - 


r 2 — 1) COS 2 ( 

72 


:= a /1 + (t 2 

— 1 ) cos : 

fdu \ 2 

1 

fdoj\ 2 

Vd~s) 

T 2 



r 2 — l) 2 cos 2 (rs) sin 2 (7 
1 + (r 2 — 1) cos 2 (rs) 


2 (r 2 — 1 ) cos (rs) sin(rs) du _ 

A 77 u. 

T OS 

On the other hand, for fixed 6 £ [0, 27t], we can apply integration by parts to deduce 
. t /t Pm r K l T 

(5.6) 


f* T du r 7T / T 

/ 2 cos(rs) sin(rs) — uds = — r / |cos 2 (ts) — sin 2 (rs)} u 2 ds. 

Jo ds J q 


Now, we use Fubini’s theorem and substitute equations (5.5) and (5.6) into (5.4). After simplifying, 
we conclude that 

1 /<9w \ 2 


l(u,u) = [ { 

J f 0 , 27 rlxr 0 , 7 r/rl ^ 


1 2 (tt) ~“ 2 }d0ds. 
[0,27r]x[0,7r/r] lT V J 


The claim follows from the change of variables x := ts and the fact that £(0, x) = ui(0, x/t). □ 


Remark 5.4. The previous lemma has an interesting consequence: the expression in polar coordi¬ 
nates of r X is a quadratic form not depending on the sphere S\(p) nor the ambient manifold M. 
This fact might suggest that all the spheres S\(p) will share the same behaviour with respect to 
stability conditions. We will illustrate this property rigorously in the main results of this section. 


Recall that any sphere S\(p) is a volume-preserving area-stationary surface, see Example 4.6. As 
we saw in Section 4.1 this implies that (A + 2AF)'(0) = 0 for any variation. Indeed, in the next 
statement we prove that S\(p) is a second order minimum of A + 2XV for certain variations which 
may preserve volume or not. 

Proposition 5.5. Let M be a 3-dimensional space form of Webster scalar curvature k. Consider 
a spherical surface S\(p) for some p £ M and A ^ 0 with A 2 + k > 0. Then, the index form of 
S\(p) satisfies l(u,u ) ^ 0 for any u £ C^iSa^)) vanishing at the poles or along the equator. As 
a consequence, S\(p) is strongly stable, in the sense that (A + 2AF) ,, (0) ^ 0 for any admissible 
variation which is C 3 off of the poles and fixes either the poles or the equator. 


Proof. Let u £ C 1 (S\(p)). From Lemma 5.3 and Fubini’s theorem, we have 


-1 /»27T ■ /»7T 

l(u,u) = - / {£f e (x ) 2 ~&(x) 2 } dx 

T Jo Ldo 


d0. 


where r := a/A 2 + k and £e(x) := y/ 1 + (r 2 — 1) cos 2 (x) u{0,x/t). Here the primes ' denote the 
derivatives with respect to x. If u vanishes at the poles, then £e(0) = £ 0 ( 77 ) = 0, and we can apply 
Wirtinger’s inequality for functions vanishing at the boundary (see for instance [40, p. 143]) to get 
Jo Ce( x ) 2 dx ^ fj fe{x) 2 dx for any 6 £ [0, 27t]. Suppose that u vanishes along the equator. Then 
u(0,tt/(2t)) = 0 for any 9 £ [0, 2ir\, and we can write t;e( x ) = cos(x)vg(x) for some continuous 
function vg : [0,7r] —> R. A straightforward computation gives us 

£' e {x) 2 - ie{x) 2 = -cos(2a :)v e (x) 2 +v' g {x) 2 cos 2 {x) - (1/2) sin(2x) (vg{x) 2 )', 
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for any x ^ 7r/2. Then, integration by parts shows that the integrals on (0, 7t/2) and (7r/2, it ) of the 
functions — (1/2) sin(2x) (vg(x) 2 )' and cos(2x) vg(x) 2 coincide. Therefore 

p2tt r pir 


X (u, u) = 


1 


v' g (x) 2 cos 2 (x) dx 


which is clearly nonnegative. The last statement follows from Theorem 5.2. 


□ 


The next result is an immediate consequence of the previous proposition. 

Corollary 5.6. A sphere S x {p) minus the poles is strongly stable, in the sense that (A + 2AF) ,, (0) 
is nonnegative for any C 3 variation. A sphere S x ( p ) removing the equator is strongly stable, in the 
sense that (A + 2XV)" (0) ^ 0 for any admissible variation which is C 3 off of the poles. In particular, 
the open hemispheres S x (p)~ and S\(p) + are strongly stable. 

Remark 5.7. Previous results related to Proposition 5.5 were obtained for the Heisenberg group 
M(0). In this particular case, it is well-known that a spherical surface S x (o) coincides with the 
union of two radial graphs defined over a horizontal disk D x and meeting along the equator. On 
the one hand, Montefalcone [42, Cor. 4.11 and App. A] showed strong stability of S\(o) for certain 
variations with velocity vector field U = (\Nh\u)N or U = uT, where it is a C 1 function with 
compact support inside an open hemisphere. On the other hand, the strong stability of S x (o) for 
admissible variations fixing the equator can be deduced from [52, Eq. (3.7)], where Ritore proved 
that S\ (o) is a global minimizer of the functional A + 2XV among finite perimeter sets H of M(0) 
such that fl C D\ x R and D x x {0} C 0. 

Remark 5.8. In spite of the previous results, the whole sphere S x (p) is not strongly stable for 
arbitrary admissible variations. To see this, it is enough to consider the deformation by Riemannian 
parallel surfaces <p s (p) ~ exp p (sN p ). This is an admissible variation by Lemma B.3 with velocity 
vector U = N. So, by taking into account Theorem 5.2, we get 

(A + 2AV)"(0)=X(1,1) = - f \N h \- 1 (l + (T 2 -l)\N h \ 2 ) 2 dS x (p)<0. 

dS\(.p) 

Clearly the variation considered in the previous example does not preserve volume. Indeed, our 
main result in this section establishes the following. 

Theorem 5.9. Let M be a 3-dimensional space form of Webster scalar curvature a. Consider a 
spherical surface S\(p) for some p £ M and X ^ 0 with X 2 + k > 0. Then, the index form of S\{p) 
satisfies X{u,u) ^ 0, for any u £ C 1 (S\(p)) with J s ^^udS\(p) = 0. As a consequence, S\{p) is 
stable under a volume constraint, in the sense that A"{ 0) ^ 0 for any volume-preserving admissible 
variation which is C 3 off of the poles. 

Remark 5.10. For a volume preserving variation of S\(p) with velocity vector U such that 
( U,N ) = u, we have f s ^^udS x (p) = 0 by (4.5) and A"{ 0) = (A + 2AV)"(0) = X(u,u) by 
(5.1). 


Proof of Theorem 5.9. By standard approximation arguments, it suffices to prove the statement 
for test functions of class C 2 . To this aim we will use double Fourier series on the rectangle 
1Z := [—7r,7r] x [—7 t/2, 7t/2]. 


Take u £ C 2 (S x (p)) with f Sx ^ udS x {p) = 0. If we define d := 6 — tt and t := x — 7t/ 2, then the 
equality in Lemma 5.3 reads 

tX(u,u) = 

where Q is the symmetric bilinear form on C l {TV) given by 




/' f dibi dib 2 
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and ip : 7?. —>■ K. is the function ip{d,t) := £($ + n,t + tt/2). Note that ip(d,— tt/2) = u(p) and 
ip{d,Tr/2) = u(p\) for any d £ [—7r,7r], whereas ip{—ir,t) = V’( 7r >*) f° r an y ^ £ [—7 t/2,7t/ 2], On 
the other hand, the mean zero condition for u together with the expression for dS\(p) in polar 
coordinates (see Lemma 3.6 (iii)) and the change of variables theorem, gives us 

0 = f £(0, x) sin(a:) dd dx = [ ip(d,t) cos(t) dd dt. 

r ‘ J [0,2-7t] x [0,7t] t ‘ JTZ 

So, in order to to prove the claim, we must show that Q(ip,ip) ^ 0 for any ip £ C 2 (1Z) such that: 


(i) ip(d, — 7t/2) and ip{d, tt/2) are constant as functions of d £ [—7r,7r], 

(ii) ip{—ir,t) = ip{ir,t), for any t £ [— 7t/2, 7t/2], 

(iii) J R ip(d,t)cos(t)dddt = 0. 


Fix ip £ C 2 (1Z) as above and define ip s (d,t) := (ip(d,t) + ip{d, — 1))/2 and ip a (d,t) := (ip(d,t) — 
ip(d,—t))/2. Clearly ip = ip s + ip a . Moreover, we have ip s (d,—t) = ip s {d,t) and ip a (d,—t) = 
— ip a (d,t ), for any t £ [— 7t/2, 7t/2] . It follows from Fubini’s theorem that Q{ip s ,ip a ) = 0 and that ip s 
satisfies equality (iii) above. In particular, Q{ip,ip) = Q(ip s ,ips) + Q(,tpa,' l Pa)- Note that 


Q(*Pa,1pa) 



ii’a)^(t) 2 } dt 


dd, 


where (ip a )ti(t) := ip a (d,t) and the derivatives are taken with respect to t. Let ( ip a )$ be the contin¬ 
uous extension of (ip a )ti to [—7r/2, 37 t/ 2] which is symmetric with respect to t = tt/2. This function 
is piecewise C 1 with mean zero. So, we can apply Wirtinger’s inequality for mean zero functions 
(see for instance [40, p. 141]) to obtain 


t /2 


I- tt/2 


(^a)tf(i ) 2 dt = 




/*37t/2 


t/2 


/ ( ' l PaY-&( t ) 2dt > OaM*) 2 * = 2 / (ipaMt ) 2 dt, 

'-tt/2 J-k/2 J — 7t/2 


and so Q{ip a ,Tpa) > 0. Hence, to prove the theorem it suffices to see that Q(ip s: ip a ) > 0. Indeed, we 
will prove that Q{ip,ip) ^ 0, for any ip £ C 2 (1Z) satisfying properties (i), (ii), (iii) and: 


(iv) ip(d, — t ) = ip{d, t), for any (d, t) £ 1Z. 


Consider a function ip as above. The Fourier series of ip with respect to the basic orthogonal 
system in L 2 (1Z) is the following sequence of functions (see [61, Sect. 7.4]): 

m n 

(5.7) S mn (d,t):= EE A ij{aij cos (id) cos(2 jt) + bij sin(*$) cos(2 jt) 

2=0 j =0 

+ Cij cos (id) sin(2 jt) + dij sin(z'd) sin(2jt)}, 

where X t j is the real constant given by 

f 1/4, for i = j = 0, 

(5.8) A ij := < 1/2, for i = 0, j > 1 or j ^ 1, j = 0, 

[l, for i.j ^ 1, 

and the coefficients aij, bij, Cij, dij are defined by 

a,ij :=E f ip(d,t) cos(id) cos(2jt) dddt, b^ := [ ip{d, t) sin(i$) cos(2jt) dd dt, 

n Jn n Jtz 

[ ip(d,t) cos(id) sm(2jt) dd dt, dij := E [ ip{d, t) sin(Z'd) sin(2jt) dd dt. 

n Jn 11 Jn 
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Note that = dij = 0 as an elementary consequence of Fubini’s theorem and (iv). On the 
other hand, the fact that ip{d,— 7r/2) = ip(d,n/2) allows us to reproduce the computations in [61, 
Sect. 3.10] to deduce that the Fourier series of dip/dt coincides with 
dS m n 

-^p(*M) = EE (—2j) Xij {a,ij cos(id) sin(2 jt) + btj sin(id) sin(2jt)}. 

2=0 j — 0 

Recall that for any cf> £ L 2 (1Z) the associated Fourier series converges to (j) in L 2 (1Z). Thus, we have 

n OO OO / n 

(5.9) / il)(d,t) h(d,t) dd dt = EE A mn \ a m n / cos(*d) cos(2 jt) h(d, t) dd dt 

J' 11 m=0 n—0 J K 


sin(zi?) cos(2 jt) h(d, t) dd dt 


m 


for any h £ L 2 (1Z). This implies the following Parseval identities for ^ and dip/dt (see [61, Sect. 7.2]): 


_2 , 
n -in 


r. OO 

/ ip 2 dddt= Y x mn(a 2 mn + b 2 mn ), 

'J 7Z ^ —n 


ra,n=0 


a /' (§.'ty 


m 


dddt= Y 4n 2 A mu{a 2 mn+b 2 mn ). 

m,n =0 

This ensures the absolute convergence of the numerical double series at the right hand sides. Now, 
the definition of the quadratic form <2(i/>, , 0) yields 


~2 2 Wb VO = E ( 4n2 “ 0 Xmn + b2 mn) = E 

m,n —0 m =0 

Here {cr m } m ^o is the sequence given by a m := a m + /3 m , where 


•— ^mO ^ m Q + ^ ^ (4n 1) Xmn ^mm •— ^mO ^mO "h ^ ^ (4?7- 1) Xmn bmn’ 

n=1 n=l 

So, to finish the proof it suffices to check that a m ^ 0 and /3 m ^ 0 for any m ^ 0. 

It is clear that /3q ^ 0 since &oo = 0. Let us see that op ^ 0. By taking into account (5.9) and 
equality j. R ip(d, t ) cos (t) dddt = 0 we get, after some easy computations 


E A 

n—0 


On 0*0n 


(~i) 


n+1 


4 n 2 — 1 


= 0, 


which is equivalent by (5.8) to 


(- 1 ) 7 

a00 — 2 2_ -—5-7 a 0n- 


An 2 — 1 


On the other hand note that, for any n £ N, 
\ 2 _ „ ^ _(—1 ) i+j a 0i a 0j 

. . E (‘ 

i,3= 


/f 2(-l)‘ O 

(Ei^nH — 4 E 


^ (4*2 - 1) (4j 2 - 1) 




E 

i,j =1 




x 0j 


(4i 2 - 1) (4j 2 - 1) (4* 2 - 1) (4j 2 - 1) 


= 4 E 


a, 


02 


= 4 (e 


X 02 


^ (4* 2 - 1) (4j 2 - 1) 

r, 2 n 2 n 

^Ez5E^E( 4 * 2 -i)' 


2=1 


4z 2 - 1 


'0i "> 


l 


E 4,-2 _ x 

7=1 J 


2=1 


where we have used that YlJLiJ 2 = ^r 2 /6. Taking limits when n —> oo and using the expression 
for aoo above we conclude that ckq ^ 0. 
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Finally, let us see that a m ^ 0 and f3 m ^ 0 for any m £ N. As ip £ C 2 (JZ) we can follow the 
arguments in [61, Sect. 3.10] to deduce that the Fourier series S mn of if in (5.7) converges uniformly 
to if in 1Z , see also [29]. Recall that if{d,'K/ 2) is a constant c € R which does not depend on 
& € [—7r, 7r]. So, we have 


OO 

E A mn {a mn (-1)” cos(wd) + b mn (-1)" sin(rm?)} = c 

m,n—0 


uniformly on [—7r, 7r]. In particular, we get also convergence in L 2 (—tt, 7r). So, for any k ^ 1, we can 
integrate on [—7r,7r] the previous equation multiplied by /i('d) := cos (kd) to obtain 


amo = 2 ^(-l) n+1 a. 

n =1 

The same computation with h(d) := sin(fc$) gives us 


& ra0 = 2^(-l)” +1 b mn 

n— 1 


The proof concludes by applying Lemma 5.11 below. 


□ 


Lemma 5.11. Let {x n } n ^o be a sequence of real numbers such that x$ = 2 1)" +1 x n and 

J])]° = i(4n 2 — 1) x„ < oo. Then, we have 



+ E (4 n 2 

n =1 


1) X n ^ 0. 


Proof. For any n £ N, we define 


-1 


n 2 n n 71 

(yy 2 (-i) i+i Xi ) + yy ^ 2 -1) x? = yy ( 4* 2 - 3 ) - 4 yy (-iy +j 


Xi Xj. 


i,j —1 

i<j 


If we show that s n Js 0 for any n G N, then the lemma follows by taking limits in the first equality 
above. Indeed, we will use induction to prove that 


(5.10) 


— 


n— 1 

E 

i=i 


n—i 2 

[(2 i - 1) Xi + 2 (-l) fc+1 %i+k + (2 n - l) 2 z 2 , n G N. 


fc=i 
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It is easy to see that (5.10) is valid for n — 1, 2. Suppose that (5.10) holds for some n ^ 2. Then 

n 

Sn+1 = S n + (4 (n + l) 2 - 3) x 2 n+1 + 4 (-l) n x n+1 Y (- 1 ) 1 x i 


i— 1 


n —1 


E [(2*-l)^ + ^ 2 (-!) 


fc+1 


%i-\-k 


+ (2n — l) 2 


k =1 


+ (4 (n + l) 2 - 3) x 2 n+1 + 4 (-1)” x n+ i E ( — 1)* 




\k +1 


2—1 
l 2 




= E[(2*-l)*i + E 2 (—1) 

2=1 fc=l 

n 

+ 4 (-l) n a; n+ i E (-1)* Xi - 4 (2n - 1) £„ x n+1 

n —1 n—i 

- Ax n+1 Y [(~l)" +i ( 2 * - 1 ) ^ + Y, 2 (- 1 ) 


+ (277- + l) 2 


°n+1 


n— 2 +/c+l 


^2+/c 


2=1 /C=l 

This gives us the desired expression for 5 n +i in (5.10) provided 

n n—1 n—i 

Y, (-1)' Xi = (-1)" (2 n - l ) x n + J2 [(-1)* ( 2 * - 1 ) Xi + Y 2 (—l) i_fe_1 x l+fe 

2=1 2=1 k— 1 

The previous equality can be checked by induction. This completes the proof. 


□ 


Remark 5.12. Montefalcone studied a 1-dimensional eigenvalues problem to prove in [42, Prop. 4.9] 
that the Pansu spheres in the Heisenberg groups HP are radially stable under a volume constraint. 
More precisely, it is shown that A" (0) Js 0 for volume-preserving admissible variations with velocity 
vector U = (|7V/j| u) N and such that u is radially symmetric. 

Remark 5.13. It is important to emphasize that our stability results in this section are valid for 
any variation of S\(p) for which formula (5.1) holds. From Theorem 5.2 this happens for admissible 
variations C 3 off of the poles. In Appendix B it is shown that the family of such variations is very 
large. However, there are other admissible variations of S\(p) with regularity less than C 3 off of the 
poles and satisfying (5.1), see Remarks A.3 for precise references. 


6. The isoperimetric problem in sub-riemannian 3-space forms 

In this section we apply our previous results to deduce some interesting consequences in relation 
to the isoperimetric problem. We first introduce some basic definitions. 

Let M be a Sasakian sub-Riemannian 3-manifold. For any Borel set Q. C M, the volume of Q. 
is the Riemannian volume V(fl) in ( M,g ). This volume functional differs by a constant from the 
signed volume defined in (4.1) since its first variation coincides, up to sign, with (4.5). As in [23] 
the perimeter of can be introduced as 

(6.1) P(fl) := sup div U dM; \U\ ^ l|, 

where U ranges over C 1 horizontal vector fields with compact support on M. Here dM and div 
are the Riemannian volume and divergence in (M, g) : respectively. If f1 is bounded by a C 2 surface 
E, then we have P(f2) = A(E) by the Riemannian divergence theorem. An isoperimetric region or 
minimizer in M is a set 12 C M satisfying P(12) ^ P(f2'), for any other 12' C M with 1/(12') = 1/(12). 
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If M is homogeneous then isoperimetric regions of any volume exist and they are all bounded 
sets, see [26, Thm. 6.1 and Lem. 4.6]. As we pointed out in the Introduction the regularity of these 
minimizers is a delicate question, whose answer is far from being established. Nevertheless, if f l is 
a C 2 isoperimetric region, then it is clear that dtd is a compact volume-preserving area-stationary 
surface. Thus, we can invoke Alexandrov’s theorem (Theorem 4.11) to get the following result, 
already known in the Heisenberg group M(0), see [54, Thm. 7.2]. 

Corollary 6.1. The boundary of any C 2 isoperimetric region f 1 in M(k) with n 0 is a spherical 
surface S\(p). 

In order to get the same conclusion in M(k) with n > 0 from Theorem 4.11, we should prove 
first that dfl is contained inside an open hemisphere of § 3 . Another possibility is to show that any 
C 2 compact volume-preserving area-stationary torus in M(«) with n > 0 (see some examples in 
Remark 4.12) is unstable under a volume constraint. In such a case, we would obtain that dfl is 
topologically a sphere by Theorem 4.4 and so, we would deduce Corollary 6.1 from Theorem 4.7. 
This is a work in progress of the authors [36]. 

The previous results might suggest that the following conjecture is probably true: in the model 
spaces M(re) any isoperimetric region is bounded by a spherical surface S\(p). In the first Heisen¬ 
berg group M(0) this conjecture was posed in 1982 by Pansu [48, p. 172] and has been supported by 
several partial results, see some references in the Introduction. Our stability result in Theorem 5.9 
provides evidence that our conjecture for any M(k) should hold. 

On the other hand, one may also ask if the spherical surfaces S\ (p) will bound isoperimetric 
regions inside any 3-space form. As happens in Riemannian geometry this is not true in general. 

Example 6.2. Consider the Sasakian flat cylinder M := M(0)/G given in Example 2.2. By con¬ 
struction, the projection map n : M(0) —> M is an isometry when restricted to the slab R 2 x (0, r). 

We compare isoperimetrically the spherical surfaces S\(p) in M with p := n(o) and A > 1/2 with 
the family of tori E# := S 1 (R) x S 1 = n(S 1 (i?) x [0, 27 t]), R > 0. Since the volume in M(0) coincides 
with the Lebesgue measure in R 3 , a straightforward computation from Lemma 3.6 (iii), see also [54, 
Rem. 7.3], implies that 

7T 2 3-7T 2 

A(S a (p)) = -, V(B x (p)) = w , 

where B\{p) is the topological 3-ball in M enclosed by S\(p). On the other hand, we have 

A( E fl ) = 4 t t 2 R, V(Q r ) = 2 t 2 R 2 , 

where flu denotes the solid 3-torus in M bounded by T,r. From the previous equalities it is easy 
to check that, for any v £ ((27/8)7r 2 ,67 t 2 ), the torus E# enclosing volume v has strictly less area 
than the sphere S\(p) of the same volume. As a consequence, none of the spheres S\(p) with 
1/2 < A < l/v^9 bounds a minimizer in M. Moreover, if we assume C 2 regularity of isoperimet¬ 
ric regions in M, then we can use the topological restriction in Theorem 4.4 together with Hopf’s 
theorem (Theorem 4.7) to deduce the existence of isoperimetric tori in M. 

In the previous example the spheres S\(p) fail to be isoperimetric when the enclosed volume 
is sufficiently large. We finish this section with the following question: are the spheres S\{p) of 
an arbitrary homogeneous Sasakian 3-manifold M boundaries of isoperimetric regions for small 
volumes? 


Appendix A. Proof of the second variation formula 

In this section we prove Theorem 5.2. We first derive a second variation formula for admissible 
variations possibly moving the singular set of arbitrary CMC surfaces in Sasakian sub-Riemannian 
3-manifolds. 
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Theorem A.l. Let E be an orientable C 2 surface, possibly with boundary, immersed in a Sasakian 
sub-Riemannian 3-manifold M. Suppose that E — Eo is C 3 and take a variation <p : I x E -A M 
which is C 3 off of Eo- Denote by U := uN + Q the associated velocity vector field, where N is the 
unit normal to E and Q := U T . IfE has constant mean curvature H and the variation is admissible, 
then the functional A + 2HV is twice differentiable at the origin, and we have 

(A.l) {A + 2HV)"(0) = [ \Nu\~ 1 {Z(u) 2 - (| B{Z) + S\ 2 + 4 (K - 1) \N h \ 2 ) u 2 } dE 

JY 

+ f &wv{{N,T){l-{B(Z),S))u 2 Z}dE 
+ f divx{(N,T)(2H\N h \u 2 -w) S}dE 

JY 

+ f divs (\Nh\W T ) dE + f divs (hiZ + h 2 S) dE, 

J Y JY 

provided all the terms are locally integrable with respect to dE. In the previous formula {Z, S'} is 
the tangent orthonormal basis defined in (2.13) and (2.14), B is the Riemannian shape operator, K 
is the Webster scalar curvature, the function w is the normal component of the acceleration vector 
field W := DjjU , and h\,h 2 are given by 

hi := 2 {H (Q, Z) + < N, T) < D S Q, Z) + |1V^| _1 (Q, S ) ((B(Z), S ) + < N, T ) 2 )} u 
+ | N h \ «Q, Z) (D S Q, S) - (Q, S > (D S Q, Z)) 

+ (N,T) < Q, Z) 2 (1 - (B(Z),S)) - (N, T) (Q, Z) (Q, S) (B(S),S), 
h 2 := -2 {H ( Q, S ) + (N, T) < D Z Q, Z) - |iV fc | (Q, Z)) u 
+ I-ZV/jI «Q, 5) (. D Z Q , Z) - (Q, Z) (D Z Q, S )) 

+ 2 H |JV„| (N, T) (Q, Z) 2 + (N, T) (Q, Z) (Q, S) (l + (B(Z),S)). 

If U = uN then the last integral in (A.l) vanishes. IfE has empty boundary and (p is supported on 
E — Eo, then all the divergence terms in (A.l) vanish. 

The theorem was proved in [55, Thm. 5.2] for some particular variations with support in E — Eo, 
velocity vector U = uN and acceleration vector W = wN. We will follow the arguments there, with 
the corresponding modifications that appear since W T and Q need not vanish on E. The idea is to 
apply differentiation under the integral sign so that, after a long calculus, we will obtain (A.l) with 
the help of Lemma A.2 below. 

Proof of Theorem A.l. We extend U along the variation by U((p s {p)) := (d/dt) |t =s <pt(p)- Let N be 
a vector field whose restriction to any E s := <^ S (E) is a unit normal vector. 

We first compute A"(0). Equation (4.2) yields 

-A(s) = / o |Jac</? s | dT,. 

is-s o 

Since the variation is admissible we can apply Lemma 2.10 to get 

(A.2) A"(0)= f {|A /l |"(0)+2|iV^| , (0)|Jac^| , (0) + |iV h ||Jac^|"(0)}dE, 

Jy-y 0 

where the derivatives are taken with respect to s. Let us compute the different terms in (A.2). 

The calculus of |Jac</3 s |'(0) and |Jacy> s |"(0) can be found in [57, Sect. 9] and [55, Lem. 5.4]. On 
the one hand, we have 

(A.3) |JacvJ s |'(0) =div s U = (-2 H R ) u + div E Q = -(2 H || + (B(S),S)) u + div E Q 
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where —2Hr := divs TV is the Riemannian mean curvature of E, and we have used that 2Hr = 
( B(Z),Z ) + ( B(S),S ) = 2H\Nh\ + {B{S),S). On the other hand, if for any p £ E — Eo, 
we choose the orthonormal basis {ei,e 2 } := {Z p ,S p } of T p E, and we take into account that 
D e U = e{u) N p — u{p) B(e) + D e Q for any e £ T p E, then we deduce 

2 

(A.4) |Jac(/? s |"(0) = divs W + (divs U) 2 + ^ |(Z? ei b r )' L | 2 

i— 1 

2 2 

-^2(R(U,ei)U,ei)- ]T (D ei U,ej) (D ej U,ei) 

i—1 i : j=l 

= div E W + |V E u| 2 + \B{Q)\ 2 + 2 (V 2 u, B(Q)> - Ric(f7, U) 

+ ( R{U , TV)E7, TV) + 2 «£> z j7, Z) <D S J7,5) - (D Z U, S) ( D S U , Z)), 

where Vem is the gradient of u relative to E and Ric is the Ricci tensor in ( M,g ). In the previous 
formula the curvature term ( R(U , N)U,N)~ Ric(U, U) can be computed from (2.3), (2.4) and (2.5). 
After simplifying, this term equals 

4(1 -K) ((TV, T ) 2 \Q\ 2 + 2 |TV h | (TV, T) (Q, S)u+ \N h \ 2 u 2 ) - \Q\ 2 - 2u 2 , 
which together with the fact that 

(D ei U,ej) = -{B(e i ),e j )u+{D ei Q,e j ), 

gives us the identity 

(A.5) |Jac^|"(0) = div E W + |V E u| 2 + |R(Q)| 2 + 2 (V 2 u, B{Q)) - \Q\ 2 

+ 2 (2 H |TV h | (B{S), S ) - (B(Z),S) 2 - l) u 2 
+ 4 (1 - K) ((TV, T ) 2 \Q\ 2 + 2 |TV h | (TV, T) (Q, S)u+ \N h \ 2 u 2 ) 

+ 2 (B(Z),S) (( D Z Q , S ) + (D S Q, Z)) u 
- 2 (B(S), S ) (D z Q, Z)u-AH |JV fc | (D S Q, S) u 
+ 2 {{DzQ, Z) ( D S Q , 5) - (DzQ, 5) (DsQ, Z)). 

Now we obtain expressions for |TV/ l | , (0) and |TV/j|"(0). From (2.16) and (2.1) it follows that 
\N h \'(s) = U(\N h \) = ( DuN , u h ) + (TV, T) (£7, Z). 

The second equality in (2.15) together with the fact that 
(A.6) D V N = -V s « - B(Q), 

implies that 

(A.7) |TV h |'(0) = —(TV,T) (5 (u) + (B(Q),S) - (Q,Z». 

Moreover, we have 

(A.8) |TV h |"(0) = ( DuDuN , + (D^TV, D v v h ) + £7((TV, T» (Q, Z) + (TV, T) £7 ((£7, Z)). 

By using (2.18), (A.6) and the third relation in (2.15), we get 
(A.9) {D V N, Duv h ) = iTV.r 1 Z{u) 2 + (TV, T) Z(u) u 

+ |TV/j| _1 ((TV, T) 2 (Q, 5) + 2 (B(Q), Z)) Z(u) 

+ |JV h | (Q, Z) S(u) + (TV, T) (B(Q), Z) u 

+ |TV /l |- 1 ((R(Q), Z) 2 + (TV, T) 2 (Q, 5) (B(Q), Z)) + |TV„| (Q, Z) (B(Q), S ). 
It is also easy to check from (2.17) that 

(A.10) U ((TV,T» = |TV h | (5(u) + (B(Q),5) - (Q,Z». 
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On the other hand 

(A.11) U((U,Z )) = (W,Z) + (D u Z,v h )(U,v h ) + (D u Z,T)(U,T) 

= (W, Z ) - (Z, DjjVh) (U, v h ) - (Z, D V T) (U, T) 

= (IV, Z) + (u + \Nh\~ 1 (N, T) (Q, S )) (Z(u) + (B(Q),Z) + (Q, S )), 

where the last equality conies from (2.18) and (2.15). It remains to compute DuDjjN. For a fixed 
point p 6 E - So we denote Ei(s) := ei(ip s ). It is clear that £^(0) = e* and [U,Ei] = 0. More¬ 
over, {Ei(s), E 2 (s)} is a basis of the tangent space to £ s at <p 3 (p). Note that {ei,e 2 , iV p } is an 
orthonormal basis of T p M. As a consequence 

2 

DuDuN = ^ (DuDuN, E t ) E t + (DuDuN, N) N. 

i= 1 

Since (AT, Ei) = 0 and [U, E,} = 0, then we get 

(DuDuN, Ei) = -2 (DuN, D v Ei) - (N, D v DuEi) = -2 £> ei [/) - (AT, DuD ei U) 

= -2 (DuN, D ei U) + <i?(J7, EOC/, AT) - (D ei W, N) 

= 2 <V S «, Z> ei J7) + 2 <B(Q), £> Bi J7) + (R(U, Ei)U, N) - (D ei W, N), 
where we have employed (A.6). Moreover, since |AT| 2 = 1 on E, we deduce 

(DuDuN, N) = -\D V N\ 2 = -|V sU | 2 - \B(Q)\ 2 - 2 (V s u,B(Q)). 

Recall that e\ = Z p and e-i = S p . Then, the previous equalities together with (S,Uh) = (N,T) lead 
us to the expression 

(A.12) < DuDuN, v h ) = (N,T) (DuDuN, S) + \N h \ (DuDuN, N) 

= (N, T) (- 2 (V E u, B{S)) u + 2 <V E u, D S Q ) - 2 (B(Q), B(S )) u 
+ 2 (B(Q),D S Q) + (R(U, S)U, N) - (D S W, N)) 

- liVftl (|Vsm| 2 + \B(Q)\ 2 + 2 <V S u, B(Q))). 

The curvature term above can be derived from (2.3), (2.4) and (2.5), so that we obtain 
(R{U, S)U, AO) = 4 (K - 1) | AO| (N, T) (Q, Z f - (Q, S ) u. 

Equations (A.12), (A.9), (A.10) and (A.11) allow us to compute |AO|"(0) from (A.8). We use the 
resulting formula together with (A.7), (A.3) and (A.5). By taking into account the equalities 

Vsu = Z(u) Z + S(u) S, B(S) = (B(Z),S ) Z + (.B(S),S ) S, 

divs Q = (D Z Q, Z) + (D S Q, S ), ( D S Q) T = (D S Q, Z) Z + (D S Q, S ) 5, 

and simplifying, we obtain 

(A.13) |AO|"(0) + 2 |AO|'(0) |Jac^r(0) + |iV h | |Jac^ s |"(0) 

= IA0T 1 Z(u) 2 + 2 (N, T) (1 - (B{Z), S )) Z(u) u + AH |IV h | (N, T) S{u) u 
+ |A0| divE W - (N, T) (D S W, N) + (N, T) (W, Z) 

+ qiu 2 + q 2 u + q 3 Z(u) + ?4 S(u) + q 5 (D Z Q, Z) + q 6 (D S Q, S) + q 7 (D S Q, Z) 

+ 2 \N h \ « D Z Q, Z) (D S Q , S) - (D Z Q, S) (D S Q, Z)) + 2 \N h \(B(Z), S ) ( D Z Q , S ) u 
+ \Nb\~ 1 (B(Q), Z) « B(Q),Z) + 2 (N, T ) 2 <Q, S)) + 2 \N h \ (Q, Z) (B(Q),S) 

+ \Nn\~ 1 (Q, S) 2 - 2 \N h \ \Q\ 2 + 4 (1 - K) \N h \ (N, T f (Q, S) 2 , 
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where the functions qi with i = 1,... ,7 are defined by 

9i = 4 H \N h \ 2 (5(5), 5 ) - 2 |JV h | (5(Z), 5) 2 + 4 (1 - K) | A h | 3 - 2 \N h \, 
q 2 = 2 (N,T) {<5(Q), Z ) (1 - (5(Z), 5)) - (Q, Z) (5(5), 5 ) 

+ 2 H |JV h | «5(Q), 5) - (Q, Z)) +4(1- A) | A,,| 2 ( Q , 5)}, 
q 3 = 2 ((A, T) (D S Q, Z) + \Nh\~ 1 ( B(Q),Z) + | A^ 1 (A, T ) 2 (Q, 5)), 
q i = 2(\N h \{Q,Z)- (A, T) (D z Q, Z», 

9s = 2 {(A, T) ((Q, Z) - (B(Q), 5)) - |JV h | <5(5), 5) «}, 

96 = 2 (A, T) (Q, Z) — 4iF |iVft,| 2 u, 

q 7 = 2 ((A, T) <5(Q), Z) + \N h \ (5(Z), 5) u). 

Note that equations (A.2) and (A.13) provide the derivative A"(0). Now, we proceed to the 
calculus of V"(0). From (4.5) we have 

V'{s)= [ (5,A)dE s = [ ((U,N)otp s ) |Jac^ s |dE, 

which implies that 

(A.14) F"(0) = /{(t/,A)'(0)+u|Jac^|'(0)}dE. 

Thus, from (A.6) and (A.3) we deduce 

(A.15) (U,N)'(0) +«|Jacy> s | , (0) = w — (Veu, Q) — (B(Q),Q) 

+ u divs Q - (2H |iV)j| + (5(5), 5)) u 2 . 

We conclude from (A.2), (A.13), (A.14) and (A.15) that 

(A.16) (A + 2HV)"(0)= f 0d£, 

JE-E 0 

where the function /3 has the expression 

(A.17) j8:= lA^- 1 Z(u) 2 + 2 (A,T) (l-(5(Z), 5}) Z(u)u 

+ 4A |iV h | <JV, T) 5(u) tt + (gi - 4A 2 |AT/,| - 2 H (5(5), 5)) w 2 
+ 2 Hw + \N h \ div E W - (A, T) ( D S W ,, AT) + (A, T) <W, Z) 

+ 92 U + (g 3 - 4FT (Q, Z)) Z(u) + (g 4 - 41F <Q, 5)) 5(it) 

+ div s (2.ffu Q) - 2 H (5(Q), Q) + g 5 <5 Z Q, Z) + g 6 <5 S Q, 5) + q 7 (D S Q, Z) 

+ 2 \N h \ ((D Z Q, Z) (D S Q, 5) - <5 Z Q, 5) (D S Q, Z)) + 2 | A h | <5(Z), 5) <5 Z Q, 5) « 

+ \Nh\~ 1 <5(Q), Z) «5(Q), Z) + 2 (A, T) 2 (Q, 5)) + 2 | A fc | (Q, Z) (5(Q), 5) 

+ |AFh,|— 1 (Q, 5) 2 - 2 |AT fc | \Q\ 2 + 4 (1 - I\) |A„| (A, T ) 2 (Q, 5) 2 . 

Let us simplify the terms containing W. From (2.16) and the fact that = (A, T) 5, it is easy to 
check that 

div E (\N h \ W T ) = |JV h | divE 1F T - <A, T) (B(S),W) + (A, T) (W, Z), 

and so 

2 Hw + \N h \ divE W - (A, T) < D S W , A*) + < A, T> (W, Z) 

= divE (|A h | W T ) - (A, T) 5(w) + (2 H (A, T) 2 - |A„| (5(5), 5)) u,. 
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Now, we can use the computations below [55, Eq. (5.20)] to infer that the first three lines of (A.17) 
equal 

(A.18) |JV fa |- 1 {^(«) 2 — (\B(Z) + S\ 2 + 4(K-l)\N h \ 2 )u 2 } 

+ div E {(TV, T) (1 - (B(Z), S )) u 2 Z} 

+ div E {(N,T) (2 H \N h \ u 2 -w)S}+ div E (|TV fc | W T ). 

To prove the statement, it suffices to show that the remainder summands in (A.17) equal divs (hi Z-\- 
I 12 S). This is a long but straightforward calculus similar to the one at the end of [55, Proof of 
Thm. 5.2]. Indeed, by using Lemma A.2 below and the formulas 

(A.19) S((N,T)) = \N h \(B(S),S), S(\N h \) = ~(N,T) (B(S),S), 

Z{(N,T » = |TV h | ((B(Z),S) - 1 ), Z(\N h \) = (N,T) (l - (B(Z),S)), 

(- D S Z , 5) = \Nh\~ 1 (TV, T) (1 + (B(Z),S)), ( D Z S , Z) = -2 H (TV, T>, 

Z{(B(Z), S )) = 4 |TV fc | (TV, T) (1 - K - H 2 ) - 2 l^" 1 (TV, T) ( B(Z),S) (l + (B(Z),S)) , 
we can deduce (A.l) from (A.16), (A.17) and (A.18). This finishes the proof. □ 

Lemma A.2. ([55, Lem. 5.5]). Let E be an orientable C 2 surface immersed in a Sasakian sub- 
Riemannian 3-manifold M. For any € C 1 (E) we have the following equalities in E — Eo 

div^cfZ) = Z(cj>) + \N h \- 1 (N,T)(l + (B(Z),S))<t>, div E (^S) = - 2H{N,T)cj>. 

Remarks A.3. 1. The regularity hypotheses in Theorem A.l are necessary to compute the different 
terms in the proof for an arbitrary variation ip. However, it is possible to derive the second variation 
for some particular variations of a surface E with less regularity and Eo = 0, see for instance [18, 
Thm. 14.5], [34, Thm. 3.7], [55, Thm. 5.2], [24, Thm. 7.3] and [27, Thm. 4.1]. In our case the C 3 
regularity of E — Eo does not suppose loss of generality since we are only interested in stability 
properties of the spherical surfaces S\(p). 

2. The fact that the variation ip is admissible is only used to differentiate under the integral 
sign in (4.2). If a variation is not admissible then A"(0) may exist or not and, in case of existence, 
it is not easy to compute it. Some examples of these situations were discussed in the Heisenberg 
group [34, Prop. 3.11] and in pseudo-Hermitian 3-manifolds, see [12, Ex. 4.3] and [24, Lem. 7.7]. 

Now, we proceed to obtain the second variation formula in Theorem 5.2. For the proof we will 
show that all the divergence terms in (A.l) vanish for variations of a spherical surface S\(p). The 
main ingredients are the integrability of |TV*,| 1 in Corollary 3.7 and a divergence theorem which 
extends the classical situation valid for vector fields supported off of the poles. 

Lemma A.4. Let X be a bounded and tangent C 1 vector field on S\(p) minus the poles such that 
divs A ( p ) X is integrable with respect to dS\{p). Then, we have 

[ div 5A(p) XdS x (p) = 0 . 

-'Sa (p) 

Proof. For any e > 0 small enough, let E £ := S\(p) — (D i(e) U H 2 (e)), where Dj(e), i = 1,2, are 
small spherical caps centered at the poles. From the divergence theorem we obtain 

2 

(A. 20 ) [ di v Sx(p) X dS x (p) = - [ ( X,rji)dl , 

2=1 ^ dDi ( e ) 

where 17 , is the inner conormal to dDi{e) in S\(p), and dl in the Riemannian element of length. Note 
that the left hand side term above tends to J s ^ ^ divg A X dS\ (p) when e —> 0 by the dominated 
convergence theorem. On the other hand, since X is bounded, the right hand side term goes to zero 
when e —> 0. This proves the claim. □ 
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Proof of Theorem 5.2. We know that S\(p) is a C°° surface off of the poles with constant mean 
curvature A. Hence, for an admissible variation <p which is also C 3 off of the poles, it is possible to 
apply Theorem A.l to compute (A + 2AF)"(0). To prove the claim it suffices to use (5.3), and to 
see that the divergence terms in (A.l) are all equal to zero. 

Consider the C 1 tangent vector field X := ( N,T ) (l — (B{Z), S)) u 2 Z defined on S\(p) minus the 
poles. Recall that (B(Z), S) = (1 — r 2 ) |iW| 2 by Lemma 3.6 (iii) and so, X is bounded. Moreover, 
by using Lemma A.2 together with some of the equalities in (A.19) and the integrability of \Nh\~ 1 
in Corollary 3.7, we deduce that divs x ( p ) X is integrable with respect to dS\(p). So, the integral of 
divs^p) X vanishes as a consequence of Lemma A.4. The same argument holds for the vector field 
X := (N, T) (2 H |1V/j| u 2 — w) S. Note that the integral of div$ A (p) (\Nh\ W T ) also vanishes by the 
divergence theorem since | N /,| W T is a Lipschitz tangent vector field on S\(p). Hence, to finish the 
proof we must show that 

(A.21) f div 5A( p) XdSxip) = 0, 

(p) 

where X := hi Z + /12 S , and the functions hi, /12 are defined below equation (A.l). At the end 
of the proof of Theorem A.l we obtained that div 5 A ( p ) X is equal to the function /3 in (A.17) by 
removing the first three lines. This long expression only contains bounded functions and the term 
1-^ft.l -1 , which is integrable on S\(p). From here, it follows that divs A (p) X is integrable with respect 
to dS\(p). Furthermore, from the definitions of hi and /12 we deduce that all the terms contained 
in X are bounded vector fields off of the poles, with the exception of 

2 liwr 1 (Q,s) « B(Z),S ) + (N,T) 2 )uZ. 

Since ( B(Z),S ) = (1 — r 2 ) lA^j 2 and \Nh\ 2 + ( N,T ) 2 = 1, then the unique unbounded vector in 
the previous formula is 2 |A^| -1 (Q,S}uZ. Hence, we cannot apply directly Lemma A.4 to infer 
equality (A.21). However, the conclusion of Lemma A.4 also holds provided the right hand side 
term in (A.20) tends to zero as e —> 0. So, to prove (A.21) we only have to see that 

(A.22) lim / |^/ i r 1 <g,5)<Z,??A U ^ = 0, * = 1,2, 

E ^-° JdDi(e) 

where Di{e), i = 1,2, are small spherical caps centered at the poles, and 7 p is a unit vector tangent 
to S\(p) and normal to dDi(e). 


Fix a positive orthonormal basis {ei, 62 } in the contact plane TL p . Let F(9, s ) := jg(s) be the flow 
of CC-geodesics of curvature A with 7 e( 0 ) = p and 7 e( 0 ) = (cos 9) ei + (sin 9) e^- As in Lemma 3.6 
we consider the unit normal along S\(p) such that Z = jg off of the poles. For e > 0 small enough, 
the curve dD i(e) is parameterized by the curve a e : [0, 27r] —!> S\(p) given by a e {9) := F{9, e). From 
equations (3.7) and (3.6), we get 

f) F 1 

a e (9) = ^-( 9,s) = V e {e) = ~(\v{e)) Z - y/v{eY + {v'{e)/2) 2 S, 

where v(e) := sin 2 (re)/T 2 and t := VA 2 + k. The previous formula implies that 

Jv{e) 2 + {v'(e)/2) 2 Z - (Xv(e)) S 

m := ra 

is a unit normal to dD i(e) tangent to S\(p). Thus, we have 

/ \N h \- 1 (Q,S)(Z, m )udl= V ^ £)2 + ^ {£)/2)2 f « Q,S)u)(a e (9))d9. 

J dDi (e) V(£) J 0 

On the one hand, it is straightforward to check that 

lim = L 
£->o v(e) 
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On the other hand, if we denote h e (9) := ((Q, S ) u)(a s (0)), then the continuity of Q and u in S\(p) 
together with equality S = ( N,T ) Vh — \ Nh \ T = — ( N,T ) J(je) — |AT/t| T , gives us \h e (9)\ < c, for 
some c > 0 not depending on e, and lim e _>o h e (0) = h(0), where 

h{9) := u(p) (( Q p , ei) sin 9 - (Q p , e 2 ) cosd). 

Thus, we can apply the dominated convergence theorem to conclude that 

/•27T /*27T 

lim / ((Q,S)u)(a e (9))d9= / h(6)d,0 = 0. 

e_>0 do do 

This yields equation (A.22) for i = 1. The case i = 2 is similar. This proves (A.21) and finishes the 
proof of the theorem. □ 


Appendix B. Examples of admissible variations 

The notion of admissible variation was introduced in Definition 5.1. For such variations we can 
apply Lemma 2.10 to differentiate under the integral sign twice in (4.2). In this appendix we will 
construct several admissible variations. We start with some immediate examples. 

Example B.l. Let U be a sub-Riemannian Killing field on M. This means that any diffeo- 
morphism <p s of the associated one-parameter group is an isometry of M. So, |Jac<y5 s | p = 1 and 
|7V/j| p (s) = |dVA|(_p) for any s € I and any p £ E. From here it is immediate that the resulting 
variation of E is admissible and the area functional is constant. Observe also that these variations 
need not fix the singular set So- 

Example B.2. Suppose that E has C 3 regular set E — Eo- Consider a variation ip which is C 3 on 
E — Eo, and such that |A^|j,(s) > 0 for any p £ E — Eo and any s £ I. Geometrically, this means 
that any ip s preserves the regular set, i.e., <p s (E — Eo) C E s — (E s )o- In particular, the function 
s i y f(s,p) := |A r ? l |p(s) |Jac<p s |p is C 2 for any p £ E — Eo, which implies conditions (i), (ii) and (iv) 
in Definition 5.1. Moreover, we have 

d 2 f 

-gjl( s ’P) = \ N h\p{s) |Jac<^ s |p + 2 \N h \' p (s) |Jac<p s | p (s) + \N h \ p (s) |Jac<p s |"(s). 

So, a possible way of constructing admissible variations is to bound the different terms in the pre¬ 
vious formula by locally integrable functions not depending on s. This is clear for |Jaccp s |p and 
|Jac<£ s |" since they are continuous on I x E. On the other hand, note that 

l^lp( s ) = ^Wp)(I N h \) = (DuN h ,v h )(ip s (p )), 

whereas 

\Nh\ p (s) = (D u D u N h ,iy h )(^p s (p)) + (DuNh, D v vh) (<£ s (p)). 

From the expression for DtjVh in (2.18) we infer that, if there is h : E —> R locally integrable with 
respect to dE, and such that \Nh\ p 1 (s) ^ h(p) for any p £ E — Eo and any s £ /, then the restriction 
of (p to /' x E is admissible for any /' CC I. In particular, if p is C 3 on E — Eo and compactly 
supported on E — Eo, then there is I 1 CC I such that the restriction of ip to /' x E is admissible. 

The analytic condition |Af/ ( |“ 1 (s) ^ h(p) : which is sufficient to get an admissible variation, was 
found by Montefalcone [43, Eqs. (25) and (26)] in the setting of Carnot groups. Since this condition 
is not easy to verify in practice, we are led to produce admissible variation in more geometric ways. 
In precise terms we will focus on normal and vertical admissible variations moving Eo- We need 
some notation and facts that will be useful in the sequel. 

Let E be an oriented C 2 surface immersed in M. For any C 1 vector field U with compact support 
on E we consider the C 1 variation ip : / x E —> M given by p s (p) := exp p (sf7 p ), where I is an 
open interval containing 0 and exp p denotes the Riemannian exponential map of M at p. Take 
a point p £ E and a vector e £ T p E. Let a : (— £q,£o) —> E be a C 1 curve with a(0) = p and 
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d(0) = e. We define the map F : (—£o,£o) x I —> M by F(e,s) := ip s (a(e)). By Remark 2.7 
we have that E(s) := (dF/de)(0,s) = e(ip s ) is a C°° vector field along the Riemannian geodesic 
7 p (s) := exp p (sU p ) satisfying [7 p ,E\ = 0 and the Jacobi equation (2.11). Note also that E( 0) = e 
and E'( 0) = D e U. Fix a basis {ei,e 2 } in T p E and let Ei(s) = ei(ip s ), i = 1,2, be the associated 
Jacobi fields. Then (i?i(s), E 2 (s)} is a basis of T 7p ( s )E and 

(B.l) |Jac <p s \ p = | E 1 x E 2 |( a ) = ({E^ 2 \E 2 \ 2 - (E 1 , E 2 )) 1/2 {s), 

which is a C°° function along j p . Moreover, the unit normal of E s := y> s (E) at <p s {p) is 


(B.2) 


N p (s) 


X E * (A 
|£i x E 2 1 


which is also C°° along 7 p . In the two previous formulas the cross product x is taken with respect 
to an orthonormal basis of the tangent space to M along 7 P (s). From the C°° regularity of |Jacc/? s | p 
and N p (s) it follows that, if |iW| p (s) > 0 for any p £ E — E 0 and any s £ /, then the conditions (i), 
(ii) and (iv) in Definition 5.1 hold. 


We are now ready to prove that the deformation of a compact surface E by means of Riemannian 
parallel surfaces is an admissible variation. 


Lemma B.3. Let E be a compact, oriented C 2 surface immersed inside a Sasakian sub-Riemannian 
3-manifold M. Then, there is an open interval f'cK such that the variation <p : I' x E —>• M defined 
by p s (p) '■= exp p (sN p ) is admissible. 


Proof. Take a point p £ E and a vector e £ T p E. Let E(s) := e(ip s ) be the associated Jacobi 
field along 7 P (s) := exp p (s7V p ). From equation (2.11) we get that (E, 7 P ) is an affine func¬ 
tion along 7 p . Thus equalities E{ 0) = e and E'{ 0) = D e N give us (E, 7 P ) = 0 along 7 p . 
As e is an arbitrary vector in T p E, the unit normal to E s at <f s {p) is Af p (s) = 7 p (s). It fol¬ 
lows that |A^| 2 (s) = 1- (j p ,T) 2 (s ) = \N h \ 2 (p) since (j p ,T) is constant along 7 P . So, we have 
|iW| p (s) |Jac| p = I^Kp) |Jac(/ 3 s | p . Hence conditions (i), (ii) and (iv) in Definition 5.1 are satis¬ 
fied. On the other hand, from equation (A.4) and the compactness of E, the derivative |Jac<p s |"(s) 
is uniformly bounded as a function of s £ I' and p £ E provided I' is a small open interval containing 
0. This finishes the proof. □ 


The previous lemma together with Example B.2 allows us to construct more general admissible 
variations based on parallel surfaces near the singular set. 

Corollary B.4. Let M be a Sasakian sub-Riemannian 3-manifold and E an oriented C 2 surface 
immersed in M with C 3 regular set E — Eo. Consider a variation ip : I x T, —> M o/E satisfying: 

(i) there is O CC E with Eo C O such that (p s (p) = exp p (sA^ p ) for s £ I and p £ O, 

(ii) the restriction of ip to (E — O) x I is of class C 3 . 

Then, there is an interval I' CC / such that the restriction of tp to I' x E is admissible. 

By a vertical variation of E we mean one of the form ip s {p) := exp p (su(p) T p ), where u £ Co(E). 
This means that the surfaces E s := y> s (51) are all vertical graphs over E. For the case u = 1 the 
variation is admissible by Example B.l since T is a sub-Riemannian Killing field and <p s coincides 
with the one-parameter group of T. In the next result we provide more examples of admissible 
vertical variations. Recall that {Z , S} is the tangent basis to E — Eo defined in (2.13) and (2.14). 

Lemma B.5. Let E be an oriented C 2 surface immersed inside a Sasakian sub-Riemannian 3- 
manifold M. Suppose that (N,T ) does not vanish along E and that |iVA, | -1 is locally integrable with 
respect to dE. Let u £ C'q(E) such that ^(u)! ^ h \Z(u)\ inYi — Eo, for some bounded function h. 
Then, the vertical variation <p s (p) '■= exp p (su(p)T p ) is admissible. 
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Proof. Take a point p £ £. The Riemannian geodesic j p (s) := exp p (s u(p) T p ) satisfies the equal¬ 
ity 7 p (s) = a p (su(p)), where a p is the integral curve of T through p. In particular, we have 
7 P (s) = zt(p) T 7 ( s ). Let {X (s), F(s)} be a positive orthonormal basis of the horizontal plane at 
7 p (s) obtained by parallel transport of a similar basis { X p ,Y p } of 7i p . For any vector e £ T p £, 
we know that the associated vector field E(s ) := e(ip s ) satisfies the Jacobi equation (2.11). By 
using (2.3) this equation reads E"(s ) + u(p) 2 E(s)h = 0. Hence, an easy integration together with 
equalities E( 0) = e and E'( 0) = e(u) T p + u(p) J(e) implies that 

(B.3) E(s) = x(s)X(s) + y(s)Y(s ) + t(s) T(s), 

where T(s) := T 7 ( s ), and the functions x(s), y{s), t(s) are given by 


(B.4) 


z(s) 

y(s) 

t(s) 


= (e,X p ) cos (u(p) s ) - (e,Yp) sin (u(p) s), 
= (e, Y p ) cos (u(j>) s ) + (e, X p ) sin (u{p) s), 
= e{u) s + (e, T p ). 


b p s 




Now suppose that p £ £ — £o and consider the orthonormal basis {ei,e 2 } of T p T, defined by 
e\ = Z p and e 2 = S p . Let Ei(s) := ei(ip s ), i = 1,2, be the associated Jacobi fields along y p . From 
equations (B.l) and (B.2), we get 

f(s,p) := \N h \ p {s) |Jac p s \ p = |(^i(s) x E 2 (s)) h \, 

where the cross product is taken with respect to the basis {JT(s), Y(s), T(s)}. Hence, a straightfor¬ 
ward computation from (B.3) and (B.4) shows that 

f(s,p) = Q p (s) 1/2 , Q p (s) := a p s 2 + i 

where 

a p := (Np,T p ) 2 Z p {u) 2 + S p {u) 2 
Note that the discriminant disc(Q p ) equals 

-41 N h \ 2 (p)(N p ,Tp) 2 Z p (u) 2 , 

which is less than or equal to 0. Moreover, the fact that (TV, T) never vanishes on £ together with 
inequality |S p (u)| ^ h{p) \Z p (u)\ implies that disc(Q p ) = 0 if and only if Z p (u ) = S p (u) = 0. From 
here we deduce that Q p (s) > 0 for any p £ £ — £o an d any s £ R. In particular, for any p € £ — £o, 
the function s i—>• f(s,p ) is C°° and so, the conditions (i), (ii) and (iv) in Definition 5.1 are satisfied. 

Finally, take p € supp(rt) D (£ — £o) and suppose a p ^ 0. Since Q p (s ) ^ Q P (—b p /(2a p )) = 
— disc(Q p )/(4a p ) and |Sp(u)| ^ h(p) \Z p (u)\ with h bounded, we obtain 


b p '.= —2 \N h \(p) S p (u), c p :=\N h \ 2 (p). 


d 2 f 


€ 


2.af 1 




disc(<5 p 

"4Q p (s)3/2 - (-disc(Q p ))V2 
((Np,T p ) 2 + h(p) 2 ) 3/2 Z p (u) 2 C((V s u) p ,Z p ) 


((N P: T P ) 2 Z p (u) 2 + S p (u) 2 ) 
\N h \(p)\(N p ,T p )\\Z p (u)\ 


3/2 


\N h \(p) 


^C'\N h \-\p), 


\N h \(p)\(N p ,T p )\ 

for some constants C, C > 0. The previous inequality also holds if a p = 0 since, in that case, b p = 0 
and ( d 2 f/ds 2 )(s,p ) = 0 for any s £ R. From the hypothesis that |TV*.| —1 is locally integrable with 
respect to d£ we conclude that condition (iii) in Definition 5.1 holds, proving the claim. □ 

Note that (N,T) ^ 0 in small neighborhoods of £o- Hence we can combine Lemma B.5 with 
Example B.2 to construct admissible variations based on suitable vertical deformations near £q. 

Corollary B.6. Let M be a Sasakian sub-Riemannian 3-manifold and £ an oriented C 2 surface 
immersed in M. Suppose that £ — £o is C 3 and the function | iV*. | 1 is locally integrable with respect 
to dY,. Consider a variation tp : I x £ —> M satisfying: 
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(i) there is a small open neighborhood O of T, o, a function u £ C 1 (0) with bounded gradient, 
and a function h bounded on O — Eo, such that 15^(it)| ^ h(p) \Z p (u)\ for any p £ O — So, 
and p s {p ) = exp (su(j>) T p ) for any s £ I and any p £ O, 

(ii) the restriction of ip to (E — O) x I is of class C 3 . 

Then, there is an interval I' CC I such that the restriction of (p to I' x E is admissible. 


The previous results provide admissible variations of a C 2 volume-preserving area-stationary sur¬ 
face E with Eo 0 inside a 3-dimensional space form. As we proved in Theorems 4.9 and 4.13, 
such a surface is either a plane C\{p) as in (3.4), a spherical surface S\(jp) as in (3.3), or a surface 
C Ml A(r) as in Example 4.14. Since all these surface are C°° off of the singular set, then Corollary B.4 
implies that C 3 perturbations of the deformation of E by Riemannian parallels are admissible vari¬ 
ations. On the other hand, the construction in Corollary B .6 only applies for C\(p) and S\(p) 
since the function || —1 is locally integrable on these surfaces, whereas it is not on C Mj ,\(r). In 
Example B.7 below we give a particular but important family of functions in S\(p) satisfying the 
condition |S(u)| ^ h \ Z(u)\, thus providing examples of admissible vertical variations. We conclude 
that the class of variations of a spherical surface S\(p) for which the second variation formula and 
the stability results in Section 5 are valid is very large. 


Example B.7. Let S\(p) be a spherical surface inside a 3-dimensional space form of curvature 
n. We consider the immersion F(9,s) := 7 g(s) introduced in Lemma 3.6, and the associated unit 
normal N. We know from Corollary 3.7 that |iW | _1 is integrable with respect to dS\(p). For any 
u £ C 1 (S\(p)) we denote u := u o F. By Lemma 3.6 (ii) and equation (3.7), we have 


z( u ) = 


S(u) = 


— 1 du , 1 1 . . du [n , 

, , ' 2 , . /0 , 9 - 'MiWKs) 

y/v(s ) 2 + (v'(s)/ 2) 2 06 OS 


where v(s) := sin 2 (ts)/t 2 and r := -\/A 2 + n. From here we deduce that, if du/dO = 0, then 
S(u) = h Z(u) for some function h bounded off of the poles. In particular, any variation ip of S\(p) 
whose restriction to small neighborhoods of the poles is of the form <p s (p) = exp p (su(p)T p ) will be 
admissible by Corollary B. 6 . Note that in a model space M(k) the hypothesis du/d9 = 0 means 
that the function u is radially symmetric, i.e., invariant under Euclidean vertical rotations. 
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